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ABSTRACT 

H^  optimal  control  theory,  based  on  singular  value  loop  sh2q)ing,  is  used  to 
synthesize  a  controller  for  the  statically  unstable  longitudinal  dynamics  of  X-29 
aircraft.  Two  design  cases  are  studied:  2-input  2-output;  and  3-input  3-output  cases. 
H„  theory  provides  a  direct,  effective  procedure  for  synthesizing  control  laws  satisfying 
specified  performance  objectives  and  robustness  specifications.  The  2  I/O  case  has 
better  performance,  a  faster  response  and  is  more  robust,  than  the  3  I/O  case. 
Discussion  and  comparison  of  results  are  given. 
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I.    INTRODUCTION 

The  main  task  for  control  system  designers  is  to  synthesize  a  control  law  which 
maintains  system  response  and  error  signals  to  within  prespecified  tolerances  while  at 
the  same  time  ensure  relatively  high  robustness.  Some  measures  of  system  robusmess 
commonly  used  are  stability  margins,  sensitivity  and  disturbance  attenuation.  TTius,  a 
robust  system  denotes  a  system  with  satisfied  system  performance,  has  large  stability 
margins,  good  disturbance  rejection  and  low  sensitivity  to  parameter  variations. 

For  the  single-input  single-output  (SISO)  systems,  some  techniques,  i.e.,  root- 
locus  method.  Bode  plot  and  Nyquist  diagram,  had  been  well  developed  to  design  and 
estimate  system  performance  and  robusmess.  Some  of  these  classical  methods  can  be 
generalized  to  multi-input  multi-output  (MIMO)  systems,  however,  not  all  of  them  do. 
For  instance,  from  a  Nyquist  diagram,  one  can  not  get  a  satisfactory  notion  of 
multivariable  stability  margins. 

In  the  past  decade,  a  lot  of  research  on  properties  of  multivariable  feedback 
system  has  been  accomplished.  Singular  value  Bode  plots  of  return  difference  and 
loop  gain  matrices  have  emerged  as  useful  indicators  of  multivariable  robusmess  [Refs. 
1,  2,  3].  Among  the  techniques  that  have  been  recently  developed,  H„,  frequency- 
weighted  LQG,  and  LQG  loop  transfer  recovery  optimal  synthesis  theories  have  made 
singular  value  loop  shaping  a  routine  matter  [Ref.  5].  For  a  multivariable  feedback 
control  system,  singular  value  loop  shaping  involves  the  manipulation  of  system  loop 
gain  over  a  specified  frequency  bandwidth  to  Lmprove  system  performance  and 
robustness. 


H„  theory  provides  a  direct,  reliable  procedure  for  synthesizing  a  controller. 
While  the  frequency-weighted  LQG  optimal  synthesis  theory  (or  "H^  theory")  and  LQG 
loop  transfer  recovery  theory  lead  to  somewhat  less  direct  but  highly  effective  iterative 
procedures  for  shaping  singular  value  Bode  plots  to  satisfy  singular  value  loop  shaping 
specification  [Ref.  5]. 

The  aim  of  this  thesis  is  to  present  a  controller  design  procedure,  bases  on  the 
H„  theory,  for  a  reduced  order,  linearized  longitudinal  dynamics  model  of  the  X-29 
aircraft  analog  backup  mode.  Some  review  about  the  properties  of  singular  value, 
sensitivity  reduction  and  return  difference  matrices  will  be  given  in  Chapter  II.  In 
Chapter  in,  H„  theory  and  weighting  matrix  are  introduced.  Application  of  H„  theory 
to  X-29  longitudinal  dynamics  model  is  made  in  Chapter  IV.  The  design  result  is 
discussed  in  Chapter  V. 


n.    FEEDBACK  PROPERTIES 

For  a  control  system,  some  properties  such  as  impulse  response  matrix,  transfer 
function  matrix  and  other  characterizations  of  the  response  to  commands  can  be 
changed  by  prefUtering  the  command  signal,  rather  than  using  feedback  loop.  But 
some  properties,  usually  termed  feedback  properties,  like  stability,  sensitivity  and 
disturbance  attenuation  can  be  altered  only  through  the  feedback  [Ref.  2].  Here,  we 
focus  on  these  feedback  properties  which  also  define  the  robustness  of  system. 

A.      RETURN  DIFFERENCE  MATRIX 

Consider  the  standard  linear  time-invariant  multivariable  feedback  control  system 

in  Fig.  2.1.    It  consists  of  the  interconnected  plant  (G),  controller  (F)  and  sensor  (H) 

forced  by  commands  (r),  measurement  noise  (n),  and  plant  disturbance  (d).     The 

prefilter  (P)  is  an  optional  element  used  to  achieve  deliberate  command  shaping. 

From  this  configuration,  we  get 

plant 

y(s)  =  G(s)u(s)  +  d(s)  (2.1) 

controller 

u(s)  =  P(s)r(s)  -  F(s)z(s)  (2.2) 

sensor 

z(s)  =  H(s)(y(s)  +  n(s))  (2.3) 

where  y(s),  d(s),  n(s),  u(s),  z(s),  r(s)  are  vectors,  and  G(s),  F(s),  H(s)  are  matrices. 

For  the  feedback  loop,  define  the  transfer  matrix 

M(s)  =  F(s)H(s) 
3 


o^ 


F 


Figure  2.1      Standard  Feedback  Control  System 


and  matrices 

L2(s)  =  G(s)M(s) 

L,(s)  =  M(s)G(s) 
are   called   return   ratio   matrices    [Ref.    2]    at   the   y(s)   node    and   the    u(s)   node, 
respectively. 
And 

I  +  L,(s) 

I  +  L,(s) 
are  return  difference  matrices  at  y(s)  node  and  u(s)  node.    Similarly,  the  y(s)  node  and 
u(s)  node  inverse  return  difference  matrices  are  defined,  respectively,  as 


I  +  L, ' 

I  +  L, ' 
Those  terms  introduced  above  are  fundamental  and  important  in  the  control  system 
synthesis.  Because  for  a  SISO  system,  the  phase  and  gain  margins  can  be  determined 
by  the  behavior  of  return  difference  as  a  function  of  frequency.  And  it  turns  out,  as 
we  will  see  later,  the  return  difference  matrix  also  provides  a  means  of  measuring 
robustness  in  a  multivariable  system. 

B.      SINGULAR  VALUES  AND  MATRIX  NORM 

In  the  design  of  a  feedback  loop,  one  needs  to  determine  the  "size"  of  matrix 
to  estimate  its  properties.  For  a  vector  X(Xi,  X2,  X3,...,xJ,  the  size  can  be  defined  as 
the  Euclidean  norm 


II  X  II  =  V  X,'  +  X,'  +  x/  +  +  x„' 

For  a  matrix,  a  more  general  concept  is  needed.     One  way  to  describe  the  size  of 
matrix  A  is  based  on  the  largest  singular  value  of  A,  as 

a(A)  =  max(  II  AX  II  /  II  X  II  ) 
where  II  •  II  denotes  Euclidean  norms.    The  singular  values  a,  of  a  rank  n  matrix  Ae 
C°"  are  defined  as  the  non-negative  square  roots  of  the  eigenvalues  of  A*A,  where  A* 
is  the  transpose  of  the  complex  conjugate  of  A.     It  is  convenient  to  order  them  as 
follow 

a,  ^  Gj  ^  Gj  ...  ^  a„ 
If  rank  r  <  n,  then  there  ait  n  -  r  zeros  singular  values,  that  is 

c^r*!  =  c^r^z  =  =  a„  =  0 


One  uses  singular  values  rather  then  eigenvalues  to  describe  matrix  size  since 
eigenvalues  can  be  a  very  poor  indicators  of  the  "size"  of  matrix  [Refs.  5,  14].  TTie 
singular  value  is  also  a  good  measure  of  the  near-singularity  of  a  matrix,  i.e.,  how 
near  the  determinant  of  a  matrix  comes  to  being  zero.  The  condition  number,  the  ratio 
of  the  largest  to  the  smallest  singular  value,  provides  information  about  sensitivity  to 
perturbation,  which  also  is  an  important  indicator  of  robustness.  Listed  as  follows  are 
some  useful  properties  of  singular  value  : 

(1)  a(A)  =  max  (II  AX  II  /  II  X  II) 

(2)  2(A)  =  min  (II  AX  II  /  II  X  II)    (the  least  singular  value  of  A) 

(3)  2(A)  <,  I  X,(A)  I  <,  a(A)  ,  where  A,  denotes  the  i-th  eigenvalue  of  A 

(4)  If  A'  exists,  2(A)  =  1  /  a(A-') 

(5)  If  A'  exists,  a(A)  =  1  /  2(A') 

(6)  a(AB)  <,  a(A)  a(B) 

(7)  Ti.,  c,'  =  Trace  (A'A)  [Ref.  5] 
In  SISO  system,  it  is  customary  to  define  sensitivity  in  terms  of  percentage 

variation  in  system  output  resulting  from  a  given  percentage  change  in  the  plant  with 
the  input  held  fixed.  A  similar  way  is  used  to  develop  the  sensitivity  function  for  a 
MIMO  system. 

C.      SENSmVITY  FUNCTION  AND  DISTURBANCE  ATTENUATION 

Consider  a  finite  dimensional  linear  time  invariant  (FDLTI)  system  in  Fig.  2.2, 
where  r,  d,  and  n  are  command,  disturbance,  and  measure  noise  resp>ectively.  If  it  is 
asymptotically  stable,  then  we  have  input-output  loop  error  and  sensitivity  relation  as 
follows  : 


Figure  2.2      Finite  Dimentional  Linear  time  Invariant  System 


input-output 

y  =  GF(I  +  GF)'(r  -  n)  +  (I  +  GF)'d  (2.4) 

loop  error 

e  =  r  -  y 
=  (I  +  GF)-'(r  -  d)  +  GF(I  +  GF)'n  (2.5) 

system  sensitivity 

AHo.  =  a  +  G'FX'AHo,  (2-6) 

where  AHo.  and  AHql  represent  the  changes  in  the  close-loop  system  and  nominal 
open  loop  system  caused  by  changes  AG  in  the  plant  G.  i.e.,  G'  =  G  -♦-  AG  [Ref.  1]. 
Here,  we  defme 

S(s)  =  a  +  L(s))-'  =  e/r  (2.7) 


R(s)  =  F(s)(I  +  L(s))-'  =  u/r  (2.8) 

T(s)  =  L(s)(I  +  L(s))'  =  I  -  S(s)  =  y/r  (2.9) 

where  L(s)  =  G(s)F(s)  is  the  loop  transfer  function,  S(s)  is  known  as  the  sensitivity 
function  matrix,  T(s)  is  complementary  sensitivity  function.  It  can  be  seen  that  loop 
transfer  function  L(s)  detemiines  the  other  three  matrices.  And,  as  will  be  seen,  the 
singular  value  bode  plots  of  these  matrices  play  important  role  in  the  design  of 
multivariable  robust  control  system. 

From  Eq.  (2.5),  it  is  seen  that  the  loop  error  "  e  "  resulting  from  disturbance 
inputs  can  be  made  small  by  making  the  sensitivity  matrix,  or  inverse  of  the  return 
difference  matrix  (I  +  GF)',  small.  From  Eq.  (2.6)  the  system  sensitivity  to  plant 
variations  may  be  also  reduced  by  decreasing  the  sensitivity  function.  One  can 
decrease  the  size  of  the  sensitivity  function  by  making 

a(S(jco))  =  a[(I  +  L(j(o))'] 
small,  or  (from  singular  value  properties  (5)) 

2(1  +  L(j(0))  (2.10) 

large. 
Note  also 

g(I  +  L(j(o))  *  g(L(jco))  for  g(L(jo)))  »  1 

This  tells  us  that  large  loop  gains  or  "  tight  "  loops  yield  good  performance.  But  one 
can  not  make  loop  gains  arbitrarily  high  over  arbitrarily  large  frequency.  Since  we  can 
see,  from  Eq.  (2.5),  large  2(L(jco))  over  a  large  frequency  range  not  only  make  errors 
due  to  r  and  d  small,  but  make  error  due  to  n  large      Large  gains  can  also  lead  to 


instability.     There  are  certain  performance  tradeoffs  between  command  /  disturbance 
error  reduction  and  sensor  noise  error  reduction  and  stability  to  be  satisfied. 

D.      STABILITY  AND  STABILFFY  MARGIN 

In  SISO  system,  gain  margin  and  phase  margin  are  the  stability  margins  that  are 
commonly  used.  Both  are  conveniently  expressed  in  terms  of  the  magnitude  of  the 
return  difference,  I  +  L(j(o).  In  MIMO  system,  the  return  difference  retains  its 
importance,  but  the  concept  of  gain  margin  and  phase  margin  become  problematic. 

The  stability  of  the  system  in  Fig.  2.2  is  directly  related  to  the  matrices  S(s), 
R(s)  and  T(s).  To  find  the  stability  margin  for  the  system  is  to  determine  a  lower 
bound  on  the  "  size  "  of  smallest  perturbation  to  the  plant  that  will  destabilize  the 
system.  And  the  singular  value  Bode  plot  of  S(s)  and  T(s)  can  be  used  as  the  measure 
of  stability  margins. 


in.    H„  CONTROL  DESIGN 

H„  control  theory  provides  a  direct,  reliable  procedure  for  synthesizing  a  feedback 
controller  designed  to  meet  singular  value  loop  shaping  requirements.  Its  capability 
includes  addressing  the  full  range  of  stability  margins,  sensitivity,  and  robustness 
optimization  and  response-shaping  problems  that  can  be  formulated  within  the  singular 
value  Bode  plot  framework.  The  standard  configuration  of  an  H„  problem  is  shown 
in  Fig.  3.1.  The  design  objective  is  to  find  a  stabilizing  controller  F(s)  for  the 
augmented  plant  P(s)  while  keeping  the  robusmess  specifications  satisfied. 
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Figure  3.1      Standard  H„  Small  Gain  Problem 
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The  purpose  of  this  chapter  is  to  present  a  brief  overview  of  H„  optimal  theory 
and  a  control  system  design  procedure  based  on  this  theory. 

A.      CONTROL  SYSTEM  FORMULATION 

As  shown  by  Eq.  (2.5),  (2.6),  and  (2.10),  the  loop  error  due  to  disturbance  and 
system  sensitivity  to  plant  variations  can  be  reduced  by  suppressing  the  sensitivity 
function  S(s).  (There  is  also  some  stability  requirements  desired  in  the  robust  system 
synthesizing.)  These  goals  can  be  achieved  by  properly  selecting  the  weighting 
functions  in  the  design  procedure. 

Shown  in  Fig.  3.2,  is  the  compensation  configuration  of  system  in  Fig.  2.2,  where 
P(s)  is  the  augmented  plant  with  weighting  functions  W,,  Wj  and  W,  which  penalizing 
the  error  signal,  control  signal  and  output  signal  respectively. 

Augmented  Plant  P(s) 
I , 

''la 
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k  w 

to  w 
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Xg-^ 

y     1 — 1 

1 
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^    1    M 

L      r 

1  ,.1 

LLT 

Ic 


y2 


Figure  3.2      Compensated  System  with  Augmented  Plant  P(s) 
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From  this  Figure,  we  get 


y2 


Jib 

Jlc 


w,   -w,G       r  u, 

0         Wj      1         u, 
0         W3G 

y^  J      Li     -g      j 

Where  the  augmented  plant  with  weighting  functions  is  denoted  as 
W,  -W,G 

0  w, 

0  W3G 

1  -G 

A  state  space  realization  of  P(s)  is  given  by 
A  B, 

P(s)  =       C.  D„ 


P(s)  = 


D, 


(3.1) 
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0 

0 
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-Do 

(3.2) 


12 


If  the  weighting  transfer  function  is  improper  (i.e.,  has  more  zeros  than  poles),  as  can 
be  the  case  for  W3,  then  no  state  space  realization  exists.    The  produce  WjG  is  in 
this  case  still  proper  and  can  be  represented  by  a  state  space  realization.     Some 
calculations  are  need  to  compute  the  state  space  realization  of  the  augmented  plant. 
From  Fig.  3.1  one  can  define  the  transfer  function  Ty,„,  as 

y.  =  Ty.^u, . 

This  transfer  function  can  be  expressed  in  terms  of  the  weighting  matrix  and  previous 
defined  loop  matrices,  i.e., 


Ui 


And  the  close  loop  transfer  function  from  input  u,  to  output  y^  is  denoted  as 


W,S 


'yi. " 

"w,s" 

y.b 

= 

W,R 

.yic . 

W3T 

'ylul 


WjR 


W3T 

Where  S,  R,  and  T  are  loop  matrices  defined  in  Eqs.  (2.7),  (2.8),(2.9). 

The  compensated  system  is  said  to  be  internally  stable  if  the  augmented  A  matrix 
of  the  compensated  system  is  stable,  i.e.,  when  external  input  signal  u,  equals  to  zero, 
the  states  of  both  P(s)  and  F(s)  will  go  asymptotically  to  zero  for  any  initial 
conditions.    And  the  controller  is  said  to  be  stabilizing. 
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B.       H„  OPTIMAL  CONTROL  PROBLEM 

The     H„  optimal  control  problem     is     defined     as     the     following  small  gain 
problem: 

For  the  given  transfer  function  matrix  P(s),  find  a  stabilizing  controller  F(s) 
such  that  the  close  loop  transfer  function  matrix  T^j^j  is  internally  stable  and  its 
infinity-norm  is  less  than  or  equal  to  1,  i.e., 

II  T^iu,  IL  <  1.  (3.9) 

Where  the  H„  -  norm  of  matrix  A  is  defined  as 

II  A  IL  =  sup  a(A(jco)).  (sup  :  the  least  upper  bound) 

This  problem  is  called  a  "  small  gain  problem  "  because  the  closed  loop  gain  is  small 
(i.e.,  less  than  or  equal  to  one).  As  mentioned  before,  the  largest  singular  values  of 
S(s)  determine  the  plant  disturbance  attenuation,  since  S(s)  is  the  same  as  the  close 
loop  transfer  function  from  disturbance  d  to  plant  output  y  (Fig.  2.2).  One  may 
specify  the  system  disturbance  attenuation  performance  with  the  frequency-dependent 
weighting  function  Wj,  that  is 

a(S(j(o))  <,  I  W.'CJco)  I 
Now,  consider  the  system  in  Fig.  3.4,  where  Aa(s)  and  A^  represent  the  additive 
plant  perturbations  and  the  multiplicative  plant  perturbations  respectively.    Defming  the 
"  size  "  of  Aa(J(0)  and  Am(J(0)  as  a(A^(j(0))  an<l  ct(_\^,(  jo)),  the  stability  fheoieins  are 
obtained  as  follows: 
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Figure  3.3      Additive  /  Multiplicative  uncertainty 

Robustness  Theorem    1: 

//  the  system  in  Fig.  3.3  is  stable  as  both  A^  and  A^,  are  zero.  Let  A^  = 
0,  the  size  of  the  smallest  stable  AJs)  for  which  the  system  become  unstable 
is 

aCAO©))  = 


1 


a(T(jco)) 
Where  T(ja))  is  the  complementary  sensitivity  function  of  the  system.  It  is  seen  that 
decreasing  the  value  of  a(T(jo)))  will  increase    the  size  of  A(j(0),  that  is,  increase  the 
stability  margin. 

Robustness  Theorem  2  : 
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Suppose  the  system  in  Fig.  3.3  is  stable  when  A^  and  A^f  are  both  zero.  Let 

A„  =  0.    Then  the  size  of  the  smallest  stable  A^(s)  which  destabilizes  the 

system  is 

1 
gCAaCjco))  =  


a(R(jo)))  [Ref.  5] 

We  can  see  that  the  smaller  a(R(jco))  is,  the  greater  will  be  the  stability  margin. 

With  properly  selected  weighting  function  matrices  W^ljco)  and  WjCjco),  the 
stability  margins  of  the  control  system  can  be  specified  by  the  following  singular  value 
inequalities, 

a(R(jco))  <  I  W;'(jco)  I 
a(T(j(o))  <  I  W3-'(j(o)  I 
Together  with  the  performance  specification  W,(jCO),  the  design  specifications  of  the 
robust  control  system  may  be  written  as 

a(S(j(o))  <  I  W,-'(j(o)  I  (3.3) 

a(R(j(o))  <  I  W,-'(j(o)  I  (3.4) 

a(T(jco))<  I  W/'CJco)  I.  (3.5) 

Thus,  for  a  robust  control  system  synthesis,  the  requirements  may  be  presented  as 
reducing  the  plant  disturbance  effect  as  much  as  possible  while  meeting  the  control 
inputs  and  the  robusmess  constraints.  Tliis  goal  can  be  achieved  by  the  singular  value 
loop  shaping  iterative  procedure  that  will  be  introduced  in  the  following  sections. 

As  a  consequence  of  singular  value  properties  and  the  definition  of  H„  norm, 
these  singular  value  inequalities  can  be  combuied  into  the  single  infinity  nonn 
specification  in  terms  of  the  close  loop  transfer  function  T  i„,  of  the  form 
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II  T,,„.  IL  ^  1 

where 

T 


^ylul 


Note,  the  matrix  Ty,„,  is  partitioned  into  submatrices  representing  performance 
requirements  and  stability  constraints. 

C.      WEIGHTING  MATRIX  SELECTION 

Since  the  plant  disturbance  can  be  attenuated  by  reducing  the  value  of  a(s(jco)) 
continuously,  a  constant  scaler  weighting  gain  y  is  added  to  W,  to  facilitate  the 
singular  value  looping  shaping  process.    That  is 

a(S(jco))  ^  \  f  *  W,-'(jo3)  I 
So,  as  we  increase  the  value  of  y,  the  sensitivity  function  is  suppressed  down 
continuously.  On  the  other  hand,  from  Eq.  (2.9)  (T(s)  +  S(s)  =  I),  we  know  that  T(s) 
goes  to  I  as  y  increases  (since  S(s)  goes  to  zero).  This  means  the  complementary 
sensitivity  is  forced  against  its  \xpptT  constraint  W3''(s)  as  the  value  of  y  increases. 
The  closed  loop  bandwidth,  therefore,  is  widened.  This  improves  the  system 
performance,  i.e.,  make  system  response  faster.  At  the  same  time,  the  singular  values 
of  closed  loop  transfer  function  Ty,„,  is  pushed  up  to  upper  limit.  The  singular  value 
loop  shaping  process  stops  as  the  condition  in  eq.  (3.9)  is  reached. 
It  is  interesting  to  note  that, 

S(s)  =  (I  +  L(s))"'  =  L(s)-' 
as  5(L(jco))  »  1,  (  i.e.,  the  system  has  large  loop  gains) 
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I.e., 

1  /  a(S(jco))  =  2(L(jco)) 
Similarly,  as  a(L(jco))  «  1, 

T(s)  =  L(I  +  L(s))-'  =  L(s) 
that  is, 

a(T(jco))  =  a(Laco)) 
These  can  be  seen  in  Fig.  3.4,  i.e.,  above  the  zero  dB  line,  2(L(j(o))  »  1,  1  /  a(S(j(o)) 
=  5(L(jco)),  while  below  the  zero  dB  line,  a(L(j(D))  «  1,  a(T(jco))  ~  a(  L(jco)).  Thus 
the  robust  system  design  consists  of  the  use  of  high  loop  gains  at  low  frequencies 
where  the  performance  specifications  are  defined  and  the  use  of  low  loop  gains  at  high 
frequencies  where  the  robustness  constraints  lie  [Refs.  5,  13].  The  system  perfomiance 
will  be  limited  within  the  range  between  a(G(s)F(s))  curve  and  g(G(s)F(s))  curve. 

The  weightings  Wj  and  Wj  are  selected  as  transfer  function  matrices  of  which 
diagonal  elements  are  frequency-dependent  constants.  And  the  size  of  these  weighting 
matrices  is  the  same  as  the  number  of  plant  output  states.  This  small  gain  problem 
is  solved  by  using  a  software  program,  i.e.,  "  hinf.m  "  provided  by  MATLAB 
ROBUST-CONTROL  TOOLBOX  package.  The  number  of  states  of  the  stabilizing 
controller  F(s)  produced  by  "  hinf.m  "  is  the  same  as  that  of  P(s).  In  addition  it  is 
required  by  "  hinf  "  that  the  D,2  submatrix  of  P(s)  (see  Eq.  (3.1))  is  full  column  rank. 
An  easy  way  to  ensure  that  this  is  the  case  is  to  choose  the  weighting  matrix  W^  with 
a  invertible  "  D-matrix  ",  e.g.,  Wj(s)  =  £l,  where  e  is  any  non-zero  number  [Ref.  5]. 
The  parameters  considered  in  selecting  these  weightings  are  gain,  order,  attenuation 
and  comer  frequency,  which  are  used  to  achieve  small  a(S(s)),  i.e.,  large  loop  gains, 
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OdB 


Figure  3.4      Singular  Value  Spec's  on  S(s)  and  T(s) 

over  the  broadest  possible  bandwidth  subject  to  the  robustness  constraints.  And  it 
should  be  noted  that  the  controller  F(s)  depends  on  the  weighting  matrices  selected. 
Once  the  weightings  are  selected,  the  control  system  design  is  in  a  computer-aided 
design  environment  and  requires  only  "  one  parameter  "  y-iteration.  The  flow  chart 
of  the  H„  iterative  procedure  is  shown  in  Fig.  3.5. 
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Figure  3.5      FL  Singular  Value  Loop  Shaping  Iterative  Process 
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IV.    H,  CONTROLLER  SYNTHESIS  FOR  X-29  FIGHTER 

The  purpose  of  this  Ch^ter  is  to  present  a  modem  aircraft  controller  design 
procedure  based  on  the  Ji.  optimal  control  theory.  The  model  used  is  the  longitudinal 
dynamic  model  of  X-29  fighter.  The  X-29  is  a  single  scat  forward  swept  wing  (FSW) 
demonstrator  aircraft  build  by  Grumman  Co.,  which  is  designed  to  be  a  new  generation 
of  tactical  fighter  with  the  advantages  of  light  weight,  low  cost  and  high  efficiency. 
The  computer-aided  design  program,  MATLAB-ROBUST  CONTROL  TOOL  BOX,  is 
used  in  controller  synthesis.  The  brief  plant  description  and  design  objective  are  given 
in  the  following  sections.  The  design  procedure  and  results  are  also  presented.  The 
script  files  specifically  written  or  modified  for  this  problem  are  listed  in  appendix  A. 

A.      PLANT  DESCRIPTION 

The  aerodynamic  advantages  of  forward  swept  wing  design  include  improved 
maneuverability,  with  virtually  spin-proof  characteristics,  better  low-speed  handling  and 
reduced  stalling  speed.  The  FSW  X-29  (see  Fig.  4.1)  was  selected  primarily  because 
its  multiple,  independently  controlled  surface  make  it  an  ideal  candidate  for  multiloop 
synthesis  of  advanced  control  mode. 

The  aircraft  is  designed  to  have  a  -35%  stability  margin  at  subsonic  speed.  A 
wing  strakc  extends  aft  fi-om  the  trailing-edge  at  each  wing  root,  each  strake  with  a 
trailing-cdge  flap  which  has  its  own  integrated  servo  actuator  to  augment  foreplanes 
for  pitch  control.  All-moving  canard  surfaces,  one  on  each  side  of  the  center  fuselage, 
outboard  of  engine  inlet  ducts  and  operated  by  servo  actuators  for  primary  pitch 
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Canard 


Figure  4.1      FSW  X-29  Fighter 

control.    In  addition  the  aircraft  is  equipped  with  flaperons  [Ref.  10]. 

The  X-29  longitudinal  dynamic  model  considered  in  this  study  is  an  analog 
reversion  mode  with  the  aircraft  trimmed  at  0.5  mark,  30,000  feet.  The  design  point 
was  one  of  the  critical  nodes  in  the  analysis  of  the  aircraft.  The  14-statc  reduced- 
order  model  includes  a  short  period  ^jproximation  of  the  aircraft  longitudinal  dynamics, 
vertical  velocity  w  and  pitch  rate  q,  and  3  fourth  order  actuator  dynamics  for  the 
longitudinal  control  surfaces  respectively,  i.e.,  the  canards,  flaperons  and  strakes. 

For  the  first  design,  a  2-input,  2-output  configuration  is  used.  Two  separated 
commands  r,  and  t^  are  input  to  the  three  control  surface  actuators  with  r,  controlling 
the  canards  and  strakes  and  Tj  controlling  the  flaps.    The  outputs  of  the  system  are  two 
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aircraft  states,  w  and  q.  The  physical  configuration  of  this  2-input,  2-output  open  loop 
actuator/aircraft  dynamics  model  is  shown  in  Fig.  4.2. 

The  second  design  example  considered  is  a  3-input,  3-output  controller  synthesis. 
Three  separated  commands  r,,  r^  and  r,  are  input  to  the  three  control  surfaces  canards, 
flaps  and  strakes  respectively.  The  output  to  be  controlled  are  w,  q  and  the  canard 
control  input  5^.  The  physical  configuration  for  this  model  is  shown  in  Fig.  4.3. 
Comparison  is  made  between  this  two  proceeding  designs. 

The  state  variables  are  listed  in  Table  4.1.  The  state  space  realization  of  the 
open  loop  plant  model  G  is  listed  in  Appendix  B.  The  poles  of  this  system  are  listed 
in  Table  4.2.  It  is  seen  that  there  is  a  positive  pole,  i.e.,  1.9550,  on  the  real  axis. 
That  means  the  X-29  has  an  unstable  short  period  mode. 

The  14-state  model  was  scaled  to  improve  the  numerical  conditioning  of  its  state 
space  representation.  The  w  state  was  transformed  to  angle  of  attack  a  using  the  initial 
forward  velocity  Uq,  i.e.,  a  =  w  /  u^.  The  units  of  the  actuator  third  derivative  states 
were  transformed  from  rad/sec^  to  le+04  rad/sec^  This  scaling  was  effective  in 
reducing  the  condition  number  of  the  system  state  matrix  A^  from  an  order  of 
magnitude    of  10'"  to  10'  [Ref.  13]. 

The  singular  value  plot  of  the  uncompensated  plant  G(s)  is  presented  in  Fig.  4.4. 
The  solid  curve  and  dashed  curve  represent  the  maximum  and  minimum  singular  values 
of  G(s)  respectively.  As  can  be  seen,  at  low  frequency,  the  open  loop  plant  has  poor 
disturbance  attenuation,  high  sensitivity  to  plant  variations  and  modeling  errors,  and  a 
small  control  bandwidth.  These  properties  will  be  improved  by  singular  vdue  loop 
shaping  procedure  presented  in  the  following  sections. 
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Figure  4.2    X-29  2-Input,  2-Output  Open  Loop  Configuration 
Source  :  W.  L.  Rogers  [Ref.  13] 
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Figure  4.3      X-29  3  I/O  Open  Loop  Configuration 
Source  :  W.  L.  Rogers  [Ref.  13] 
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Table  4.1    Discription  of  The  Uncompensated  X-29  state  variables 


State        Description  Units 

a  angle-of-attack  rad 

q  pitch  rate  rad/sec 

S^  canard  control  input  rad 

6^  flap  control  input  rad 

<5g  strake  control  input  rad 

6^  canard  control  rate  rad/sec 

6^  flap  control  rate  rad/sec 

^g  strake  control  rate  rad/sec 

6^  canard  control  accel.  rad/sec 

6^  flap  control  accel.  rad/sec^ 

S^  strake  control  accel.  rad/sec^ 

6^  canard  control  jerk  le+04  rad/sec^ 

'S\  flap  control  jerk  le+04  rad/sec^ 

S'^  strake  control  jerk  le+04  rad/sec 

Source  :  W.  L.  Rogers  [Ref.  13] 
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Table  4.2    Uncompensated  X-29  open  loop  poles 


poleg  = 


■2.2 
•2.2 

•1.4 

•1.4 

1.9 


-1 

•2 

■5, 

-5, 

•5, 

•5, 

•5, 

■2, 

•2 


746e+02 
746e+02 
455e+02 
491e+02 
550e+00 


031e+02 
155e+00 
518e+01 
518e+01 
506e+01 
506e+01 
067e+01 
172e+01 
115e+01 


3201e-f02i 
3201e+02i 


8255e+01i 
8255e+01i 
8410e+01i 
8410e+01i 
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Figure  4.4      Singular  Value  Plot  of  The  Uncompensated  Plant  G(s) 
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B.      DESIGN  SPECIFICATIONS  AND  WEIGHTING  FUNCTIONS 

Since  the  design  objective  is  to  attenuate  the  disturbance  effect  as  much  as 

possible      within  loop   gain  constraint   and  robustness   specifications.      The  design 

specifications,  in  this  case,  are  : 

Performance  Spec: 

Minimize  the  sensitivity  function  as  much  as  possible. 

Robustness  Spec: 

1).     Attenuate       the    close-loop    singular    values    of   the    complementary 

sensitivity  matrix  by  20  db  at  frequencies  beyond  co  =  100  rad/sec 

(15.9  Hz)  (This  will  ensure  that  the  system  will  have  sufficient  stability 

margin  to  tolerate  variations  the  loop  transfer  function  of  magnitude  as 

large  as  a  factor  of  ten  at  frequency  (O  =  1(X)  rad/sec). 

2).     a  second-order  roll-off  beyond  CO  =   100  rad/sec  (This  ensures  the 

controller  is  proper  and  matches  the  uncontrolled  drop-off). 

With  these  considerations,  we  have  covered  the  most  important  issue  in  feedback 

control  system  design.     Once  the  design  specifications  are  decided,  the  selection  of 

weighting  matrices  is  the  next  step.    For  the  2-input  2-output  system,  the  weighting 

matrices  are  chose  as  foUows  : 

.01  (  1  +  s  /  .01  ) 

(YW,(s))-'  =  y'  ♦      ♦  I  (4.1) 

(  1  +  s  /  100  )  (2x2) 


W2(s)  =  -  0.025  *  I.  (4.2) 

(4  X  4) 
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1000 

W;\s)  =   *  I 

s'  (2  X  2) 


(4.3) 


to  penalize  the  error  "  e  ",  control  input  "  u  "  and  output  "  y  "  respectively. 

The  singular  value  plot  of  W,"'(s),  W/'Cs)  weighting  functions  is  given  in  Fig. 
4.5.    The  robustness  specification  W/'Cs)  has  20  dB  drop-off  and  40  dB  roll-off  at 
(0  =  100  rad/sec  which  satisfied  the  design  requirements.    As  mentioned,  W,'(s)  and 
W,''(s)  should  be  so  chose  to  ensure  the  0  dB  crossover  frequency  of  Wj  is  sufficiently 
lower  than  the  0  dB  crossover  frequency  of  Wj'Cs).    WjCs)  is  selected  to  penalize  the 
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control  input  u,  and  the  weighting  factor  0.025  is  chosen  to  confine  the  actuators' 
response  within  their  design  limits.  The  selected  Wils)  also  ensures  that  the  submatrix 
D,2  of  P(s)  has  full  column  rank  as  required  in  the  two  Riccati  equation  solution  to  the 
H„  problem.  The  H„  theory  requires  that  D,2(s)  matrix  is  full  column  rank  in  order 
for  an  optimal  H„  controller  to  eliminate  infinite  energy  inputs.  If  we  decrease  the 
weighting  factor,  more  energy  will  input  to  the  system,  in  this  case,  the  actuators' 
response  may  exceed  the  limits  [Ref.  13]. 

C.      DESIGN  PROCEDURE 

Before  starting  the  singular  value  loop  shaping  procedure,  it  is  necessary  to  find 
the  balanced  state  space  realization  for  the  unaugmented  plant  for  a  better  numerical 
condition.  This  can  be  seen  from  the  different  condition  numbers  of  state  matrix  Ag 
before  and  after  balanced  realization,  which  are  5.526e+04  and  2.5567e+02  respectively. 

1.      2-INPUT,  2-OUTPUT  DESIGN  CASE 

With  the  selected  weighting  matrices,  the  augmented  plant  P(s)  of  the 
balanced  plant  is  of  order  eighteenth,  W,(s)  and  W^  each  adding  two  states  to  the 
system  G(s).  Since  WjCs)  is  not  a  proper  transfer  function  Eq.  (4.3)  (i.e.,  it  has  more 
zeros  than  poles),  so  it  has  no  state  space  realization.  But  W3(s)G(s)  is  proper  and  can 
be  realized  in  state  space  form. 

With  the  computer-aided  design  program,  y  is  chosen  as  1  for  the  first  try 
of  the  iteration  process.  As  continuously  increasing  the  value  of  y,  the  singular  value 
of  cost  function  Ty,„,  is  pushed  up  to  its  upper  limit  gradually  until  no  solution  exists 
for  any  larger  y.  After  several  iterations,  y  equiiJs  to  2.6,  and  the  maxiinum  singular 
values  of  Ty,„,  reached  the  all  pass  limit  (i.e.,  0  dB  line).    Fig.  4.6  is  the  singular  value 
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Figure  4.6      Cost  Function  Ty,„,  for  y  =  2.6 

plot  of  cost  function  Tj,,„,,  the  solid  line  represents  the  maximum  singular  values  and 
dashed  line  represents  the  minimum  singular  values.  Fig.  4.7  and  4.8  arc  the  singular 
value  plots  of  complementary  sensitivity  and  sensitivity,  the  dashed  line  and  dotted  line 
again  represent  the  maximum  and  minimum  singular  values  respectively.  The 
complementary  function  is  pushed  flat  against  the  specified  limit  W,',  as  the  value 
of  gamma  is  increased,  the  closed  loop  bandwidth  is  widened  (to  30  rad/sec)  and  the 
response  speed  is  increased.  Similarly,  the  sensitivity  function  is  pushed  down  against 
its  limit  y'W,'  as  y  value  increased.  The  singular  value  of  I  +  G(s)F(s)  loop  gains 
are  increased  to  about  50  dB  at  low  frequency  for  this  design,  which  markedly 
improves  the  X-29  performance  properties.  This  large  loop  gains  indicate  good 
disturbance  attenuation  and  low  sensitivity  to  plant  variations  over  a  control  bandwidth 
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Figure  4.7      Complementary  Sensitivity  Function  T(s)  and  W,"'  for  y  =  2.6 
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of  about  3  rad/sec  and  satisfy  the  stability  requirements.  Note,  in  Fig.  4.8,  W,'(s)  = 
a(S(jo)))  at  low  frequency.  And  there  is  about  40  dB  drop-off  from  (O  =  10  rad/sec  to 
CO  =  0.1  rad/sec  which  is  the  range  of  interesting.  The  state  space  realization  of  the 
18-state  stabilizing  controller  F(s)  is  listed  in  Appendix  C.  Two  uncontrollable  states 
are  removed  using  PRO-MATLAB  function  "  minreal.m  ".  A  fast  transient  poles  is 
eliminated  while  retaining  the  most  important  characteristics  of  the  closed  loop  system. 
Thus,  the  final  controller  has  15  states.  The  poles  of  this  reduced-order  controller  are 
listed  in  Table  4.3. 


Table  4.3    Poles  of  The  Reduced-Order  Controller 


poleacpr  = 

-2. 1954e-l-02 

-6.5170e-(-01    -I-  9.5444e+01i 

-6.5170e-i-01    -  9.5444e-i-01i 

-1.4487e-l-02 

-1.3577e-t-02 

-5.3024e-i-01    +  4  .  9271e-i-01i 

-5.30246-1-01    -  4.9271e-i-01i 

-4.1572e-l-01    -I-  3  .  1587e-(-01i 

-4.1572e+01  -  3.1587e+01i 

-1.31806-1-01  -I-  1.9602e-l-01i 

-1.31806-1-01  -  1.96026-1-011 

-3.01186-1-01  +  1.83616-l-Oli 

-3.01186-t-Ol  -  1.836l6-»-01i 

-1.00006-02 

-1.00006-02 
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Figure  4.9      X-29  Feedback  Configuration 

Fig.  4.9  shows  the  closed  loop  system  of  compensated  X-29,  where  F(s)  is 
the  15th  order  controller  and  G(s)  is  the  14th  order  X-29  plant.  Since  F(s)  is  placed 
in  series  with  plant  G(s),  the  command  vector  r  is  composed  of  the  reference 
commands  r,  and  r^  to  controlled  outputs  a  and  q.  So  ,  the  closed  loop  compensated 
system  has  2  inputs,  2  outputs  and  29  states.  The  poles  of  this  reduced-order  closed 
loop  system  are  listed  in  Table  4.4.  The  singular  value  plot  of  output,  inverse-return 
difference  matrix  I  -t-  (G(s)F(s))"'  (also  called  the  output,  multiplicative  return  difference 
matrix)  for  the  compensated  system  is  shown  in  Fig.  4.10.  The  minimum  singular 
value  of  this  matrix,  i.e.,  <T„i„[I  +  (G(jco)F(jco))'l.  provides  a  measure  of  the  aircraft's 
gain  and  phase  margins  with  respect  to  multiplicative  modeling  errors  using  the 
universal  gain  and  phase  margin  curve  (Fig.4.11)  [Ref.  4]. 
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Table  4.4    Poles  of  The  2  I/O  Closed  Loop  System 


-2.2739e+02  +  2.3194e+02i 

~2.2739e+02  -  2.3194e+02i 

-2.1723e+02 

-6.4876e+01  +  9.6310e+01i 

-6.4876e+01  -  9.6310e+01i 

-1.3558e+02 

-1.4595e+02 

-1.4469e+02  +  9.2844e-02i 

-1.4469e+02  -  9.2844e-02i 

-1.1118e+02 

-4.9782e+01  +  4.9587e+01i 

-4.9782e+01  -  4.9587e+01i 

-5.2188e+01  +  4.8038e+0li 

-5.2188e+01  -  4.8038e+01i 

-5.2464e+01  +  4.8363e+01i 

-5.2464e+01  -  4.8363e+01i 

-4.1968e+01  +  3.2086e+01i 

-4.1968e+01  -  3.2086e+01i 

-4.1045e+01 

-2.7932e+01  +  1.3447e+01i 

-2.7932e+01  -  1.3447e+01i 

-1.1994e+01  +  1.1959e+01i 

-1.1994e+01  -  1.1959e+01i 

-1.9603e+00 

-2.6697e+00 

-3.8935e+00 

-1.0123e+01 

-1.6232e+01 

-1.9727e+01 


It  is  seen  in  Fig.  4.10  that  a^„[I  +  (G(s)F(s))']  drops  to  about  -3  dB  at 
frequencies  between  0.3  rad/sec  and  20  rad/sec.  Entering  the  vertical  axis  of  the 
universal  curve  at  singular  value  of  0.71  (i.e.,  -3  dB).  The  gain  and  phase  margins  of 
compensated  X-29  near  0  dB  crossover  frequency  are  -11  dB  to  +5  dB  and  ±42  deg 
respectively.  This  is  more  stable  than  the  -8  dB  to  +4  dB  and  ±35  deg  gain  and  phase 
margins  for  the  typical  design  of  a  fighter  aircraft. 
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Figure  4.10    Singular  Value  Plot  of  The    Output  Inverse-Return  Difference 
Matrix  for  2  I/O  System 


t^.ind  +  G'^JW)) 


^.in(I+Q(JW)) 


PHASE  MARGIN  ±8 


10  0  IS 

GAIN ^MARGIN    (dB) 


20 


Figure  4.11     Universal  Gain  and  Phase  Margin  Curve 
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Shown  in  Fig.  4.12  are  the  closed  loop  Bode  plots  of  a  and  q  respectively. 
Tlie  solid  curves  are  the  response  to  r,  and  dashed  curves  are  to  r^.  It  is  seen  that,  at 
low  frequency,  a  is  "all  pass"  (i.e.,  has  0  dB  gain)  to  input  1  but  has  very  small 
response  (-50  dB)  to  input  2.  Similarly,  as  seen  that  q  is  "all  pass"  to  i^  but  has  small 
response  to  r,.  This  indicates  the  decoupling  of  a  from  q  to  r,  and  i^.  Note,  the 
system  performance  at  high  frequencies  satisfy  the  robustness  requirements  for  both 
cases. 

The  time  response  of  a  and  q  to  1  deg/sec  (0.01745  rad/sec)  step  input  are 
presented  in  Fig.  4.13  and  4.14.  It  is  seen  from  Fig.  4.13  that  the  response  of  a  is  fast 
while  with  rise  time  of  0.3  seconds  and  0.1  second  to  r,  and  r^  resj)ectively.  We  also 
can  see  that  the  decoupling  of  a  and  q  to  the  two  commands  from  both  plots.  If  we 
reduce  the  weighting  function  W2(s),  more  energy  will  be  input  to  system,  and  that 
energy  increases  the  response  speed  and  the  decoupling  will  be  more  pronounced. 

Referring  to  Fig.  4.15  and  4.16,  the  time  response  of  control  deflections  5^, 
bf,  and  5..  It  is  seen  that  the  peak  magnitudes  are  0.013  rad,  0.006  rad,  and  0.007  rad 
for  5^,  5f,  and  5,  respectively,  which  are  well  within  the  X-29's  control  surface 
deflection  limits  listed  below  : 

(1).  canards  (leading  edge)  :  30  deg  up  /  60  deg  down, 

(2).  flaps  (trailing  edge)  :  10  deg  up  /  25  deg  down, 

(3).  strakes  (trailing  edge)  :  30  deg  up  and  down. 
The  peak,  control  rates  (see  Fig.  4.17  and  4.18)  are  0.0013  rad/sec,  0.0004  rad/sec,  and 
0.0008  rad/sec  respectively,  which  are  much  smaller  than  the  X-29  actuator  minimum 
design  requirements  of 
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X-29  alpha  RESPONSE  TO  0.01745  rad  /  1   sec  STEP  INPUT 
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Figure  4.13    Step  Response  of  a  for  2  I/O  System 
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Figure  4.14    Step  Response  of  q  for  2  I/O  System 
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Figure  4.15    Step  Response  of  Canard  Control  Input  5^  and  Flap  Control 
Input  5f 
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Figure  4.16    Step  Response  of  Strake  Control  Input  5, 
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xlO-*  X-29  DFDOT  FOR  0.01745  rad  /  1  sec  STEP  INPUT 
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Figure  4.18    Step  Response  of  5f  and  6, 
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(1).  canards  :  1.75  rad/sec 

(2).  flaps  :  0.87  rad/sec 

(3).  strakes  :  0.52  rad/sec. 
The  response  is  well  within  the  actuator  limits  but  one  also  has  a  very  small  step- 
input. 

2.      3-INPUT,  3-OUTPUT  DESIGN  CASE 

The  second  design  exanq)le  we  use  is  the  same  X-29  dynamic  model  as  in 
the  previous  example,  but  with  different  input  and  output  vectors,  i.e.,  r,,  r^,  and  t^  to 
control  canards,  flaps  and  strakes  separately  and  outputs  are  a,  q,  and  5^  (see  Fig. 
4.3).    Thus,  it  is  a  3-input,  3-output  system  with  14  states. 

The  same  design  specifications  are  used  for  this  case,  and  the  weighting 

functions  are  chose  as  follows  : 

.01  (  1  +  s  /  .01  ) 

(YW,(s))-'  =  y'  *    '  *  I  (4.4) 

(  1  +  s  /  100  )  (3x3) 


W,(s)  =  -0.025  *  I.  (4.5) 

(6  X  6) 


1000 

W3-'(s)  =    *  I  (4.6) 

s'  (3  X  3) 


The  augmented  plant  P(s)  is  an  20th  order  system  with  "W|(s)  and  W;(s)  each  adding 
three  states  to  the  plant  G(s).  Following  the  same  iterative  procedure  as  of  the 
previous  example,  the  cost  function  reaches  its  all  pass  limit  as  y  =  1.3  (see  Fig.  4.19), 
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COST  FUNCTION  Tylul  (Gamma  =  1.3) 
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Figure  4.19    Cost  Function  Ty,u,  for  y  =  1.3 

then  the  iterative  process  stops.  The  solid,  dashed,  and  dotted  curves  in  this  plot 
represent  the  maximum,  medium,  and  minimum  singular  values  of  Ty,u,  respectively. 
The  controller  F(s)  produced  has  20  states,  that  is  the  same  as  number  of  states  of 
P(s).  Three  states  arc  removed  using  "  minreal.m  "  function  and  one  state  is 
eliminated  by  "  obalreal.m  "  program.  So,  the  final  controller  has  16  states,  the 
controller  /  plant  series  has  30  states  in  total  (uncompensated  plant  G(s)  has  14  states). 
The  poles  of  this  reduced-order  closed  loop  system  are  listed  in  Appendix  E. 
Comparing  Fig.  4.20  with  Fig.  4.7  and  4.8  shows  that  this  3  I/O  system  are  not  as 
robust  as  the  2  I/O  model.  The  3  I/O  system  has  smaller  disturbance  attenuation, 
larger  sensitivity  to  plant  variations  and  modeling  errors,  smaller  control  bandwidth  and 
closed  loop  bandwidth.    Fig.  4.21  is  the  singular  value  plot  of  output,  inverse-return 
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X-29  SV  PLOT  (I  +  inv(GF)) 


160 


FREQUENCY  -  rad/sec 

Figure  4.21    Output,  Inverse  Return  Difference  Matrix  for  3  I/O  System 

difference  matrix.  It  is  seen  that  the  c^„[  I  +  ((G(s)F(s))'  ]  (i.e.,  the  dotted  curve) 
drops  to  approximately  -6  dB  between  0.2  rad/sec  and  13  rad/sec.  Enter  the  vertical 
axis  of  the  universal  curve  (Fig.  4.11)  at  singular  value  of  0.5  (i.e.,  -6  dB),  the  gain 
and  phase  margins  near  0  dB  crossover  frequency  are  -6  dB  to  +3  db  and  ±  28  deg 
respectively.  It  indicates  that  the  stability  margins  of  the  3  I/O  system,  with  respect 
to  multiplicative  modeling  errors,  are  less  than  the  values  desired  of  a  fighter  aircraft, 
i.e.,  -8  dB  to  +4  dB,  ±35  deg  gain  and  phase  margin  respectively.  Obviously,  it  is 
also  not  as  stable    as  the  2  I/O  system  discussed  before. 

From  Fig.  4.22,  the  Bode  plots  of  a  and  q  to  the  3  inputs,  we  also  can  see 
that  a  is  decoupled  from  q  to  r,  and  r,.    And  it  is  interesting  to  see  that  both  a  and 


q  have  very  Little  response  to  Vy    This  also  can  be  seen  from  Fig.  4.23,  i.e.,  r. 


has 
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Figure  4.23    Step  Response  of  a  and  q  for  3  I/O  System 
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almost  no  effect  on  both  a  and  q.    This  indicates  that  there  is  not  much  difference 
to  system  response  between  2  inputs  design  or  3  inputs  design. 

Shown  in  Fig.  4.24  and  4.25,  the  peak  values  of  control  deflections  6^,  6,, 
and  6,  (i.e.,  0.005  rad,  0.0017  rad,  and  0.015  rad  respectively)  are  all  well  within  the 
design  limits  listed  before,  and  are  also  smaller  than  those  of  the  2  I/O  system  since 
the  loop  gains  are  smaller.  Similar  results  are  obtained  from  the  plots  of  control  rates 
(Fig.  4.26  and  4.27),  i.e.,  all  responses  stay  within  their  design  limits. 

Thus,  we  may  conclude  that  the  performance  of  this  3  I/O  control  system  satisfies 
the  design  specifications  but  is  not  better  than  those  of  2  I/O  system  designed 
previously. 


49 


XlO-3 


9 

c 

■a 

s 
C 
s 


-4 


Response  of  DC  to  0.01745  rad/sec  Step  Input 


0.5 


TO  r3 


TO  rl 


1.5  2  2.5 

TIME  -  sec 


3.5 


xlO-3 


Response  of  DF  to  0.01745  rad/sec  Step  Input 


1.5 


-2 


to  rl 


0.5 


1.5 


2.5 


Time  -  ■ec 
Figure  4.24    Step  Response  of  6,  and  5f  for  3  I/O  System 
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Figure  4.25  Step  Response  of  6,  for  3  I/O  System 
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Figure  4.27  Step  Response  of  6,  and  5.  for  3  I/O  System 
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V.    CONCLUSIONS  AND  RECOMMENDATIONS 

Although  the  mathematic  work  behind  the  H„  theory  is  very  complicated,  the 
implementation  of  this  method  is  easy  using  a  well  developed  computer-aided  design 
software.  From  the  design  procedure  presented  above,  it  is  seen  that  H„  method 
provides  a  direct,  systematic  and  effective  way  to  synthesize  controllers  for  the 
statically  unstable  longitudinal  dynamics  of  X-29  aircraft.  The  two-command  control 
design  example  demonstrates  a  fast  response,  large  stability  margins,  good  disturbance 
attenuation  and  low  sensitivity  performance.  While  the  3-input  control  configuration 
shows  a  somewhat  less  performance  due  to  the  smaller  loop  gains. 

It  is  seen  in  both  design  cases  that  the  two  control  outputs,  angle-of-attack  a 
and  pitch  rate  q,  are  decoupled  to  the  command  inputs.  This  is  an  important  feature 
which  is  needed  to  effect  the  advanced  control  modes  in  which  aircraft  flight  path  and 
pitch  attitude  are  independendy  controlled.  Since,  as  shown  in  the  3-input  case,  the 
separated  commanded  strake  has  almost  no  influence  on  the  control  outputs,  we  may 
conclude  that  a  reasonable  and  effective  controller  can  be  based  on  the  use  of  two- 
command  control  in  the  X-29  control  system  . 

In  both  designs,  we  found  the  peak  values  of  control  surface  deflections  for  the 
actuators  and  control  rates  are  all  much  smaller  than  their  design  limits.  The  selected 
W2  weighting  matrix  ensures  this  compatibility  but  at  a  severe  penalty  in  performance. 

Further  reduced  order  modeling  of  the  controller  is  possible  but  was  not 
accomplished  in  this  thesis. 


53 


APPENDIX  A 
Matlab  Script  Files 


diary  hx29 

dat 

format  short  e 

dispC 

') 

dispC 

This  script  file  is  designed  to  solve  the  Hinf  optimal  control  ') 

dispC 

problem  for  the  X-29.   The 

14th  order  FDLTI  model,  in  state  ') 

dispC 

space  form,  is  that  of  the 

X-29  aircraft  and  actuator  dynamics.') 

dispC 

Two  states  are  those  of  the  aircraft  dynamics,  i.e.,  alpha  ') 

dispC 

and  q.   The  remaining  12  states  are  the  dynamics  of  the  ') 

dispC 

three,  fourth  order  actuators,  i.e.,  the  canard,  flaperon,  and  ') 

dispC 

strake  actuators.   The  order  of  the  unbalanced  states  is  as  ') 

dispC 

follows: ') 

dispC 

') 

dispC 

alpha,  q,  dc,  df,  ds,  dcdot,  dfdot,  dsdot,  dcdbldot,') 

disp( ' 

dfdbldot,  dsdbldot,  dctrpldot,  dftrpldot,  dstrpldot') 

dispC 

') 

disp( ' 

Given  the  open  loop  transfer  function  G (s) =Cinv (Is-A) B+D,  a  ') 

dispC 

stabilizing  controller  F(s) 

will  be  found  such  that  the  Hinf  norm' 

dispC 

of  Tylul  is  less  than  or  equal  to  one.   ') 

dispC 

') 

dispC 

') 

%pause 

clc 

dispC 

X-29  aircraft  and  actuator 

state  space  representation') 

dispC 

') 

dispC 

') 

ag=[-.4181d+00   .9960d+00   -.2269d-01 

-.1213d+00   -.1948d-01   -.9493d-03  ... 

.4427d- 

04 

-.6712d-04   .1451d-05   -.2162d-04   -.3540d-05   0.0   0.0   0.0; 

.5474d+01   -.3424d+00    .2585d+01 

-.1386d+01   -.1058d+01   .3898d-02  ... 

-.1164d 

-01   -.6397d-02   -.2509d-03   - 

■.5362d-03   -.2912d-03   0.0  0.0  0.0; 

0.0 

0 

0   0.0   0.0   0.0   1.0   0.0 

0.0   0.0   0.0   0.0   0.0   0.0   0.0; 

0.0 

0 

0   0.0   0.0   0.0   0.0   1.0 

0.0   0.0   0.0   0.0   0.0   0.0   0.0; 

0.0 

0. 

0   0.0   0.0   0.0   0.0   0.0 

1.0   0.0   0.0   0.0   0.0   0.0   0.0; 

0.0 

0 

0   0.0   0.0   0.0   0.0   0.0 

0.0   1.0   0.0   0.0   0.0   0.0   0.0; 

0.0 

0. 

0   0.0   0.0   0.0   0.0   0.0 

0.0   0.0   1.0   0.0   0.0   0.0   0.0; 

0.0 

0. 

0   0.0   0.0   0.0   0.0   0.0 

0.0   0.0   0.0   1.0   0.0   0.0   0.0; 

0.0 

0. 

0   0.0   0.0   0.0   0.0   0.0 

0.0   0.0   0.0   0.0   l.Oe+04   0.0   0.0; 

0.0 

0 

0   0.0   0.0   0.0   0.0   0.0 

0.0   0.0   0.0   0.0   0.0   l.Oe+04   0.0; 

0.0 

0 

0   0.0   0.0   0.0   0.0   0.0 

0.0   0.0   0.0   0.0   0.0   0.0   l.Oe+04; 

0.0 

0 

0   -.1479d+04   0.0   0.0   -. 

1143d+03   0.0   0.0   -.2529d+01   0.0  ... 

0.0   - 

.2697d+03   0.0   0.0; 

0.0 

0 

0   0.0   -.1491d+04   0.0   0. 

0   -.1149d+03   0.0   0.0   -.2536d+01  ... 

0.0   0 

.0 

-.2701d+03   0.0; 

0.0 

0 

0   0.0   0.0   -.5302d+05   0. 

0   0.0   -.1816d+04   0.0   0.0  ... 

-.1790d+02   0.0   0.0   -.6053d+03] 

bg=[0.0 

0 

0 

0.0 

0 

0 

0.0 

0 

0 

0.0 

0 

0 

0.0 

0 

0 

0.0 

0 

0 

0.0 

0 

0 

0.0 

0 

0 

0.0 

0 

0 

0.0 

0 

0 

0.0 

0 

0 

.1479d+04  0.0; 
0.0  .1491d+04; 
.5308d+05   0.0] 
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cg= [1.0   0.0 
0.0   1.0 


0.0 
0.0 


0.0 
0.0 


0.0 
0.0 


0.0 
0.0 


0.0 
0.0 


0.0 
0.0 


0.0 
0.0 


0.0 
0.0 


0.0 
0.0 


0.0 
0.0 


0.0 
0.0 


0.0; 

0.0] 


dg=zeros (2 ) 

Ipause 

dispC     ') 

disp('     Balanced  realization  of  the  X-29  state  space  representation') 

dispC     ') 

[agbl ,bgbl,cgbl ,g, t]=obalreal (ag, bg,cg) 

ag=agbl;  bg=bgbl;  cg=cgbl; 

%pause 

w=logspace(-3 , 3 , 100)  ; 

svg=sigTna  (ag,  bg,  eg,  dg,  1 ,  w)  ;       svg=2  0*logl0  (svg)  ; 

semilogx (w, svg) 

title('X-29  SV  PLOT  OF  PLANT  G') 

xlabelC  FREQUENCY  -  rad/sec') 

ylabel ( ' SV  -  dB' ) 

grid 

meta  hx29gsv 


disp 
disp 
disp( 
disp( 


poleg=eig(ag) ,   tzerGg=tzero(ag, bg, cg,dg) 


disp( 
disp( 
disp( 
disp( 
disp( 


condag=cond (ag) ,   rcondag=rcond (ag) 


disp( 

disp( 

%pause 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp  ( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 


') 

Calculate  the  poles  and  transmission  zeros  of  the  balanced  open' 

loop  plant' ) 
') 


') 
') 
Determine  determine  the  condition  number  of  ag') 

') 
') 


') 
') 


') 


') 


1) 
') 


<<  Design  Specifications  >>   ') 

Robustness  Spec.  :  -40  dB  roll-off,  -20  db  g  100  Rad/Sec.') 
Associated  Weighting:') 


') 


2) 


') 


-1        1000  ') 

W3(s)  =    *  I        (fixd)') 

2      2x2') 
s        ') 


Performance  Spec:  minimizing  the  sensitivity  function') 

as  much  as  possible.') 
Associated  Weighting:') 


') 
-1        -1  .01(100s  +  1) 
Wl(s)  =  Gam   * * 

(.is  +  1) 


') 
I') 

2x2') 
') 


') 


where  "Gam"  in  this  design  is  iteratively  updated  from  1') 
w=logspace (-3,3,100) ; 
k=1000;   mn=[2  2];  tau=0.0; 
nuw3i  =  [0.0  k];  dnw3i  =  [1.0  0  0]; 
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svw3i  =  bode(nuw3i,dnw3i, w) ;  svw3i  =  20*logl0 (svw3i) ; 

nuwli  =  [1.0  0.01];  dnwli  =[0.1  1.0]; 

svwli   =   bode(nuwli, dnwli, w) ;    svwli   =   20*loglO (svwli) ; 

aw2=-0.025*eye(2) ;    bw2=zeros (2 ) ;    cw2=2eros (2 ) ;    dw2=-0. 025*eye (2) ; 

dispC       ') 

dispC       ') 

disp('  (strike  a  key  to  see  the  plot  of  the  weightings  ...)') 

Ipause 

semilogx(w, svwli, w, svw3i) 

grid 

title ('X-29  Design  Specifications') 

xlabel ( 'Frequency  -  rad/sec') 

ylabel('l/Wl  &  1/W3  -  db') 

text(.01,0, '1/Wl(s) ') 

text( .5, 100, '1/W3 (s) ') 

neta  hx29specs 

%pause 

clc 

disp 

dispC  ') 

disp 

disp 

disp('  1) 

dispC   ') 

disp('  2) 

dispC   ') 

disp 

disp 

disp 

disp 

disp 

dispC  ') 

disp('      where  ') 

disp 

disp('  Tylul  = 

disp 

disp 

dispC   ') 

dispC   ') 

dispC   ') 

disp 

Ipause 

clc 

dispC   ') 

dispC  ') 

disp 

dispC   ') 

disp 

disp 

disp 

disp 

disp 

disp 

disp 

disp 

disp 

disp 

disp 

disp 

dispC   ') 


<<  Problem  Formulation  >>') 

Form  an  augmented  plant  P(s)  with  these  two  weighting  functions:') 

Gam*Wl  penalizing  error  signal  "e"') 

W3  penalizing  plant  output  "y"' ) 

and  find  a  stabilizing  controller  F(s)  such  that  the  Hinf  norm') 
of  TF  Tylul  is  less  than  or  equal  to  one,  i.e.') 

|Tylul|    <  or  =  1, ') 
F(s)        inf  ) 


-1 

') 

Gam*Wl*(I    +   GF) 

= 

-1 

W3*GF*(I    +   GF) 

') 

Gam  *  Wl  *  S 
W3  *  (I  -  S) 


(strike  a  key  to  continue 


') 
') 


)') 


<<  DESIGN  PROCEDURE  >>') 

***************************  *') 

[Step  1].  Do  plant  augmentation  (run  AUGMENT. M  or    *') 

AUGX29.M) 


Balance  the  augmented  plant  for  better 
numerical  condition  if  necessary 


[Step  2] 

[Step  3].  Do  Hinf  synthesis  with  "Gam"  =  1 


[Step  4].  Redo  the  plant  augmentation  for  a 
higher  "Gam"  and  rerun  HINF.M 


*') 
*') 
*') 
*') 
*') 
*') 
*') 
*') 
*') 


*********** 


*************' 
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dispC   ') 

disp('  (strike  a  key  to  continue  ...)') 

%pause 

clc 

dispC   ') 

dispC   ') 

dispC   ') 

disp('  Assign  the  cost  coefficients  "Gam"  with  Gani=l  ') 

dispC       ') 

disp('  serving  as  the  baseline  design  ....') 

dispC   ') 

disp( '  ' ) 

disp('  augment    %  Plant  augmentation  of  the  X-29  dynamics') 

disp( '  ' ) 

dispC    ') 

dispC    ') 

haugx29 

dispC   ') 

dispC   ') 

disp('  (strike  a  key  to  continue  ...)') 

%pause 

clc 

dispC    ') 

dispC    ') 

%disp('      Do  state  space  balancing  on  the  augmented  plant  if  needed') 

%disp('    ') 

%disp('    ') 

%[abal,bbal,cbal,g,t]=obalreal(A, [Bl  B2 ] , [CI ;C2 ] ) 

%A=abal,    Bl=bba^ ( : , 1 : 2) ,    B2=bbal ( : , 3 : 4 ) ,    Cl=cbal ( 1 : 4 , : ) ,    C2=cbal (5: 6 , : ) 

dispC         ') 

dispC         ') 

disp('    The  transmission  zeros,  poles  and  condition  number  of  the  augmented') 

disp('     plant  follow.   In  addition,  determine  if  (A,B1)  &  (A,B2)  are  ') 

disp('     stabilizable  and  if  (CI, A)  &  (C2,A)  are  detectable.') 

dispC    ') 

dispC    ') 

t2eroaug=tzero(A, [Bl  B2 ] , [CI ;C2 ] , [Dll  D12;D21  D22]),  poleaugA=eig (A) 

condaugA=cond (A) ,  rcondaugA=rcond (A) 

eps=eps 

toldef=10*max(size (A) ) *norm(A, 1) *eps 

tol=100*eps*norm( [A  Bl]) 

[Alc,Blc,Clc,t,k]=ctrbf (A,Bl,Cl,tol) 

tol=100*eps*norm( [A  B2 ] ) 

[A2c,B2c,C2c,t,k]=ctrbf (A,B2,C2,tol) 

tol=100*eps*norm( [A;C1] ) 

[Alo,Blo,Clo,t,k]=obsvf (A,Bl,Cl,tol) 

tol=100*eps*norTn(  [A;C2] ) 

[A2o,B2o,C2o,t,k]=obsvf (A,B2,C2,tol) 

clear  condag  rcondag  poleg  tzerog  svw3i  svwli 

clear  condaugA  rcondaugA 

clear  functions 

%pause 

dispC     ') 

dispC     ') 

disp( '    ' ) 

disp('     hinf     %  Running  script  file  HINF.M  for  Hinf  optimization') 

dispC    ') 

aretype= ' Schur ' 

hinf 

dispC   ') 
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dispC 

') 

dispC 

clear  : 

functions 

Ipause 

dispC 

') 

dispC 

') 

dispC 

State 

dispC 

(acp, 

disp( ' 

') 

(strike  a  key  to  continue  ...)') 


space  representation  of  the  full  order  controller') 
bcp,  ccp,  dcp)  with  its  poles  and  condition  number') 

acp , bcp , ccp , dcp 

polec=eig (acp) 

condacp=cond (acp) 

dispC     ') 

disp('    Minimal  realization  of  the  controller') 

dispC     ') 

toldef=10*max(size (acp) ) *norm(acp, 1) *eps 

tol=100*eps*norm( [acp   bcp;ccp  dcp]) 

[acpm, bcpm, ccpm, dcpm]=minreal (acp, bcp, ccp,  dcp) 

dispC  ') 

disp('     Balanced  realization  &  model  reduction  of  the  minimal  controller' 

dispC     ') 

[acpbl , bcpbl , ccpbl , g, t ]=obalreal (acpm, bcpm, ccpm) 

elim=[10] 

[acpr , bcpr , ccpr , dcpr ]=modred (acpbl , bcpbl , ccpbl , dcpm, elim) 

dispC     ') 

disp('     Poles,  controllability,  observability,  and  condition  of  the  ') 

disp('     balanced,  reduced  order  controller') 

dispC     ') 

poleacpr=eig (acpr) 

tol=100*eps*norm ( [acpr  bcpr]) 

[acpc, bcpc,ccpc, t,k]=ctrbf (acpr , bcpr , ccpr , to 1) 

tol=100*eps*norm( [acpr ; ccpr] ) 

[acpo, bcpo, ccpo, t , k]=obsvf (acpr, bcpr , ccpr, tol) 

condacpr=cond (acpr) ,   rcondacpr=rcond (acpr) 

acp=acpr;  bcp=bcpr;  ccp=ccpr;  dcp=dcpr; 

dispC     ') 

disp('     CLTF  Tylul  (acl,  bcl,  ccl,  del)  and  its  poles  (reduced  order)') 

[acl, bcl , ccl , dcl]=lf tf (sysp, dimp, acp, bcp, ccp, dcp) 

polet=eig (acl) 

%pause 

dispC    ') 

dispC    ') 

hpltopmod  %  Preparing  singular  values  for  plotting 

end 

dispC    ') 

dispC    ') 

disp('  (strike  a  key  to  continue  ...)') 

%pause 

dispC    ') 

dispC    ') 

disp('    Open  loop  state  space  representation  of  controller/plant  series') 

dispC    ') 

[algf , blgf ,clgf ,dlgf ]=series(acp, bcp, ccp, dcp, ag, bg,cg,dg) 

polol=eig (algf) 

dispC    ') 

dispC    ') 

disp('  (strike  a  key  to  continue  ...)') 

%pause 

dispC    ') 

dispC    ') 
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diary  h3x29.dat 

format  s 

hort  e 

dispC 

h3x29.ni') 

dispC 

') 

dispC 

This  script  file  is  designed  to  solve  the  Hinf  optimal  control  ') 

dispC 

problem  for  the  X-29.   The  14th  order  FDLTI  model,  in  state  ') 

dispC 

space  form,  is  that  of  the  X-29  aircraft  and  actuator  dynamics.') 

dispC 

Two  states  are  those  of  the  aircraft  dynamics,  i.e.,  alpha  ') 

dispC 

and  q.   The  remaining  12  states  are  the  dynamics  of  the  ') 

dispC 

three,  fourth  order  actuators,  i.e.,  the  canard,  flaperon,  and  ') 

dispC 

strake  actuators.   The  order  of  the  unbalanced  states  is  as  ') 

dispC 

follows: ' ) 

dispC 

') 

dispC 

alpha,  q,  dc,  df,  ds,  dcdot,  dfdot,  dsdot,  dcdbldot,') 

dispC 

dfdbldot,  dsdbldot,  dctrpldot,  dftrpldot,  dstrpldot') 

dispC 

') 

dispC 

Given  the  open  loop  transfer  function  G (s) =Cinv (Is-A) B+D,  a  ') 

dispC 

stabilizing  controller  F(s)  will  be  found  such  that  the  Hinf  norm' 

dispC 

of  Tylul  is  less  than  or  equal  to  one.   ') 

disp( ' 

') 

dispC 

') 

%pause 

clc 

dispC 

X-29  aircraft  and  actuator  state  space  representation') 

dispC 

') 

dispC 

') 

ag=[-.4181d+00   .9960d+00   -.2269d-01   -.1213d+00   -.1948d-01   -.9493d-03  ... 

.4427d- 

04   -.6712d-04   .1451d-05   -.2162d-04   -.3540d-05   0.0   0.0   0.0; 

.5474d+01   -.3424d+00    .2585d+01   -.1386d+01   -.1058d+01   .3898d-02  ... 

-.1164d 

-01   -.6397d-02   -.2509d-03   -.5362d-03   -.2912d-03   0.0  0.0  0.0; 

0.0 

0.0   0.0   0.0   0.0   1.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0; 

0.0 

0.0   0.0   0.0   0.0   0.0   1.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0; 

0.0 

0.0   0.0   0.0   0.0   0.0   0.0   1.0   0.0   0.0   0.0   0.0   0.0   0.0; 

0.0 

0.0   0.0   0.0   0.0   0.0   0.0   0.0   1.0   0.0   0.0   0.0   0.0   0.0; 

0.0 

0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   1.0   0.0   0.0   0.0   0.0; 

0.0 

0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   1.0   0.0   0.0   0.0; 

0.0 

0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   l.Oe+04   0.0   0.0; 

0.0 

0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   l.Oe+04   0.0; 

0.0 

0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   0.0   l.Oe+04; 

0.0 

0.0   -.1479d+04   0.0   0.0   -.1143d+03   0.0   0.0   -.2529d+01   0.0  ... 

0.0   - 

.2697d+03   0.0   0.0; 

0.0 

0.0   0.0   -.1491d+04   0.0   0.0   -.1149d+03   0.0   0.0   -.2536d+01  ... 

0.0   0 

.0   -.2701d+03   0.0; 

0.0 

0.0   0.0   0.0   -.5302d+05   0.0   0.0   -.1816d+04   0.0   0.0  ... 

-.1790d 

+02   0.0   0.0   -.6053d+03] 

bg= [0.0 

0.0   0.0 

0.0 

0.0   0.0 

0.0 

0.0   0.0 

0.0 

0.0   0.0 

0.0 

0.0   0.0 

0.0 

0.0   0.0 

0.0 

0.0   0.0 

0.0 

0.0   0.0 

0.0 

0.0   0.0 

0.0 

0.0   0.0 

0.0 

0.0   0.0 

.1479d-t-04   0.0   0.0; 

0.0 

.1491d+04   0.0; 

0.0 

0.0   .53( 

)8d+05] ; 
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cg=[1.0  0.0 
0.0  1.0 
0.0   0.0 


0.0 

0.0 
1.0 


0.0 
0.0 
0.0 


0.0 
0.0 
0.0 


0.0 
0.0 
0.0 


0.0 
0.0 
0.0 


0.0 
0.0 
0.0 


0.0 
0.0 
0.0 


0.0 

0.0 
0.0 


0.0 
0.0 
0.0 


0.0 
0.0 
0.0 


0.0  0.0; 
0.0  0.0; 
0.0   0.0]; 


dg=zeros (3,3)  ; 

Ipause 

dispC     ') 

disp('     Balanced  realization  of  the  X-29  state  space  representation') 

dispC     ') 

[agbl , bgbl,cgbl,g, t]=obalreal (ag, bg,cg) 

ag=agbl,  bg=bgbl,  cg=cgbl, 

w=logspace (-3,3, 100) ; 

svg=sigTOa (ag, bg, cg,dg, 1, w) ;   svg=2  0*logl0 (svg)  ; 

seinilogx(w,  svg) 

title ('X-29  SV  PLOT  OF  PLANT  G') 

xlabelC  FREQUENCY  -  rad/seC) 

ylabel('SV  -  dB') 

grid 

%meta  h3x29gsv 


disp( 
disp( 
disp  ( 
disp( 


poleg=eig(ag) 


disp( 
disp( 
disp( 
disp  ( 
disp( 


condag=cond(ag) ,   rcondag=rcond (ag) 


disp( 

disp( 

%pause 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp( 

disp  ( 

disp( 

disp( 

disp( 


'} 

Calculate  the  poles  and  transmission  zeros  of  the  balanced  open') 
loop  plant' ) 
') 


tzerog=tzero (ag, bg, eg, dg) 


') 
') 
Determine  determine  the  condition  number  of  ag' 

') 
') 


') 
') 


') 


1) 


') 


<<  Design  Specifications  >>   ') 

Robustness  Spec.  :  -40  dB  roll-off,  -20  db  g  100  Rad/Sec.') 
Associated  Weighting:') 


-1 
W3(s)  = 


1000  ') 

*  I        (fixd) ') 

2      3x3') 
s        ') 


') 


2) .  Performance  Spec. 


') 


minimizing  the  sensitivity  function') 
as  much  as  possible.') 
Associated  Weighting:') 


') 
-1        -1  .01(100s  +  1) 
Wl(s)  =  Gam   * * 

(.is  +  1) 


') 
I') 

3x3') 

') 


where  "Gam"  in  this  design  is  iteratively  updated  from  1') 


w=logspace(-3, 3 , 100) ; 


60 


k=1000;   inn=[3  3];  tau=0.0; 

nuw3i  =  [0  0  k];  dnw3i  =  [l.O  0  0]; 

svw3i  =  bode (nuw3i,dnw3i, w) ;  svw3i  =  20*logl0 (svw3i) ; 

nuwli  =  [1.0  0.01];  dnwli  =[0.1  1.0]; 

svwli  =  bode(nuwli,dnwli,w) ;  svwli  =  20*loglO (svwli) ; 

aw2=-0.025*eye(3) ;  bw2=2eros (3) ;  cw2=zeros (3) ;  dw2=-0 . 025*eye(3) ; 

dispC   ') 

dispC   ') 

disp('  (strike  a  key  to  see  the  plot  of  the  weightings 

%pause 

semilogx (w, svwli , w, svw3i) 

grid 

title('X-29  Design  Specifications') 

xlabel ( 'Frequency  -  rad/sec') 

ylabel('l/Wl  &  l/W3  -  db') 

text(.01,0, '1/Wl(s) ' ) 

text(.5, 100, '1/W3(s) ') 

%meta  h3x29spc 

%pause 

clc 

disp 

disp 

disp 

disp 

disp 

disp('   ') 

disp 

disp('   ') 

disp 

disp('   ') 

disp 

disp 

disp 

disp 

disp 

disp 

disp 

disp('  (strike  a  key  to  continue  ...)') 

%pause 

clc 

disp 

disp 

disp 

disp 

disp 

disp 

disp 

disp 

disp 

disp 

disp 

disp 

disp 

disp 

disp 

dispC   ') 

dispC   ') 

disp 

%pause 

clc 


<<  Problem  Formulation  >>') 

Forin  an  augmented  plant  P(s)  with  these  two  weighting  functions:') 

1).  Gam*Wl  penalizing  error  signal  "e"') 

2).  W2  penalizing  control  input  "u"' ) 

2).  W3  penalizing  plant  output  "y"') 

and  find  a  stabilizing  controller  F(s)  such  that  the  Hinf  norm') 
of  TF  Tylul  is  less  than  or  equal  to  one,  i.e.') 

|Tylul I    <  or  =  1, ' ) 
F(s)        inf) 


<<  DESIGN  PROCEDURE  >>') 


[Step  1].  Do  plant  augmentation  (run  H3AUGX29.M  ) 

[Step  2].  Balance  the  augmented  plant  for  better 
numerical  condition  if  necessary 

[Step  3].  Do  Hinf  synthesis  with  "Gam"  =  1 

[Step  4].  Redo  the  plant  augmentation  for  a 
higher  "Gam"  and  rerun  HINF.M 


*') 
*') 
*') 
*') 
*') 
*') 
*') 
*') 
*') 
*') 
*') 


(strike  a  key  to  continue 


)') 
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dispC   ') 

dispC   ') 

dispC   ') 

disp('  Assign  the  cost  coefficients  "Gam"  with  Gain=l  ') 

dispC       ') 

disp('  serving  as  the  baseline  design  ..,.') 

dispC   ') 

disp( '  ' ) 

disp('  augment     %  Plant  augmentation  of  the  X-29  dynamics') 

disp( '  ' ) 

dispC    ') 

dispC    ') 

h3augx29 

dispC   ') 

dispC   ') 

disp('  (strike  a  key  to  continue  ...)') 

Ipause 

clc 

dispC    ') 

dispC    ') 

%disp('     Do  state  space  balancing  on  the  augmented  plant  if  needed') 

IdispC    ') 

%disp('    ') 

%[abal,bbal,cbal,g,t]=obalreal(A, [Bl   B2] , [C1;C2] ) 

%A=abal,    Bl=bbal ( : , 1 : 2) ,    B2=bbal ( : , 3 : 4) ,    Cl=cbal (1 : 4 , : ) ,    C2=cbal (5 : 6 , : ) 

dispC         ') 

dispC         ') 

disp('    The  transmission  zeros,  poles  and  condition  number  of  the  augmented'; 

disp('    plant  follow.   In  addition,  determine  if  (A,B1)  &  (A,B2)  are  ') 

dispC     stabilizable  and  if  (CI, A)  &  (C2,A)  are  detectable.') 

dispC    ') 

dispC    ') 

tzeroaug=tzero(A, [Bl  B2 ] , [CI ;C2 ] , [Dll  D12;D21  D22]),  poleaugA=eig(A) 

condaugA=cond (A) ,  rcondaugA=rcond (A) 

eps=eps 

toldef=10*max(size (A) ) *norm(A, 1) *eps 

tol=100*eps*norm( [A  Bl]) 

[Alc,Blc,Clc,t,k]=ctrbf (A,Bl,Cl,tol) 

tol=100*eps*norm( [A  B2]) 

[A2c,B2c,C2c,t,k]=ctrbf (A,B2,C2,tol) 

tol=100*eps*norm( [A;C1] ) 

[Alo,Blo,Clo,t,k]=obsvf (A,Bl,Cl,tol) 

tol=100*eps*norm( [A;C2] ) 

[A2o,B2o,C2o,t,k]=obsvf (A,B2,C2,tol) 

clear  condag  rcondag  poleg  tzerog  svw3i  svwli 

clear  condaugA  rcondaugA 

clear  functions 

%pause 

dispC     ') 

dispC     ') 

disp( '    ' ) 

disp('      hinf     %  Running  script  file  HINF.M  for  Hinf  optimization') 

disp( '     ' ) 

aretype= ' Schur ' 

hinf 

dispC   ') 

dispC   ') 

disp('  (strike  a  key  to  continue  ...)') 

clear  functions 

Ipause 
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dispC 
dispC 
dispC 
dispC 
dispC 


') 


State  space  representation  of  the  f 
(acp,  bcp,  ccp,  dcp)  with  its  poles 


') 


ull  order  controller') 
and  condition  number') 


acp , bcp , ccp , dcp 

polec=eig(acp) 

condacp=cond (acp) 

dispC     ') 

disp('     Minimal  realization  of  the  controller' 

dispC     ') 

toldef=10*max (size (acp) ) *norm(acp, 1) *eps 

tol=100*eps*norm( [acp  bcp;ccp  dcp]) 

[acpm, bcpm, ccpm,dcpm]=minreal (acp, bcp, ccp, dcp 

dispC  ') 

disp('     Balanced  realization  &  model  reduct 

dispC     ') 

[acpbl , bcpbl, ccpbl,g, t]=obalreal (acpm, bcpm, cc 

elim=[10] 

[acpr , bcpr, ccpr, dcpr]=modred (acpbl , bcpbl , ccpb 

dispC     ') 

disp('     Poles,  controllability,  observabili 

disp('    balanced,  reduced  order  controller 

dispC     ') 

poleacpr=eig (acpr) 

tol=100*eps*norm( [acpr  bcpr]) 

[acpc, bcpc, ccpc, t , k]=ctrbf (acpr, bcpr, ccpr, tol 

tol=100*eps*norm( [acpr;ccpr] ) 

[acpo, bcpo, ccpo, t , k]=obsvf (acpr, bcpr , ccpr , tol 

condacpr=cond (acpr) ,   rcondacpr=rcond (acpr) 

acp=acpr;  bcp=bcpr;  ccp=ccpr;  dcp=dcpr; 

dispC     ') 

dispC     CLTF  Tylul  (acl,  bcl,  ccl,  del)  and 

[acl , bcl, ccl,dcl]=lftf (sysp,dimp, acp, bcp, ccp, 

polet=eig(acl) 

%pause 

dispC    ') 

dispC    ') 

h3plot 

end 

dispC    ') 

dispC    ') 

dispC 

%pause 

dispC    ') 

dispC    ') 

dispC 

dispC 

[algf , blgf , clgf ,dlgf ]=series (acp, bcp, ccp, dcp, 

polol=eig(algf ) 

dispC 

dispC 

dispC 

%pause 

disp( ' 

dispC 

dispC 

dispC 

dispC 

dispC    ') 


>) 

ion  of  the  minimal  controller') 

pm) 

1 ,dcpm,elim) 

ty,  and  condition  of  the  ') 
) 


its  poles  (reduced  order)') 

dcp) 


%  Preparing  singular  values  for  plotting 

(strike  a  key  to  continue  ...)') 

Open  loop  state  space  representation  of  controller/plant  series') 

ag,bg,cg,dg) 

(strike  a  key  to  continue  ...)') 


) 

) 

Closed  loop  state  space  representat 

controllability,  observability,  and 

loop  acgf  matrix.  ') 


ion  of  controller/plant  series,') 
condition  number  of  the  closed') 
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dispC    ') 

[acgf , bcgf , ccgf , dcgf ]=f eedbk (algf,blgf,clgf,dlgf,2) 

tol=100*eps*nonn  ( [acgf  bcgf]) 

[acgf c, beg f c, ccgf c,t,k]=ctrbf (acgf, bcgf , ccgf , tol) 

tol=100*eps*norTn(  [acgf , -ccgf  ] ) 

[acgf o, bcgf o, ccgf o,t,k]=obsvf (acgf, bcgf, ccgf , tol) 

condacgf=cond (acgf ) 


dispC 

') 

dispC 

') 

dispC 

(St 

%pause 

dispC 

') 

dispC 

Poles  of  the  closed  loop  system' 

disp( ' 

') 

polcl=ei 

■  g(acgf) 

h3aly 

%h3rsp 

end 

(strike  a  key  to  continue  ...)') 
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%  Plant  Augmentation  for  the  X-29  Hinf  problem  as  W3  is  not  a 
\   proper  transfer  function.   Includes  contingency  for  adding  W2  to 
%  ensure  dl2  is  full  column  rank.   This  script  file  is  designed  for 
%  the  X-29  system  without  theta  as  a  state,  ie,  3  inputs  &  3  outputs. 
%disp('    ') 


') 


%disp( ' 
%disp(' 

%Gam=gama ( 1 , i) 

Gam  =  input ( ' 

cgb  =  l/k*[cg(l, 

cg(2, 

eg  (3, 

dgb  =  l/k*[cg(l, 

eg  (2, 

eg  (3, 


<<  Plant  Augmentation  >>') 


•Gam"  =  ') 


*ag*ag; 
*ag*ag] ; 


Input  the  cost  coefficient 

: ) *ag*ag; 

: ) *ag*ag*ag*tau+cg (2 , : 

: ) *ag*ag*ag*tau+cg ( 3 , : 

: )*ag*bg; 

: ) *ag*ag*bg*tau; 

: ) *ag*ag*bg*tau] ; 
nwl  =  Gam*[dnwli;0  0;0  0;0  0;dnwli;0  0 ; 0  0;0  0;dnwli]; 
dwl  =  nuwli; 

sysw2=[aw2  bw2;cw2  dw2 ] ;  xw2=3; 

%[A,B1,B2,C1,C2,D11,D12,D21,D22 ]=haugmod (ag, bg, eg, cgb, dg,dgb, nwl ,dwl ,mn) 
[A,B1, B2 ,C1,C2 ,D11,D12 , D21 , D22 ] =haugx29pl (ag, bg, eg, cgb, dg, dgb, nwl , dwl , sysw2,xw2, 


tdispC 

IdispC 

%disp(' 

% 

% 


') 


-  -  State-Space  (A, B1,B2 , CI ,C2 ,D11 ,D12 , D21,D22)  is  ready  for') 
the  Small-Gain  problem  -  -  -') 


End  of  AUGX29.M 
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disp('     Closed  loop  state  space  representation  of  controller/plant  series,') 

disp('     controllability,  observability,  and  condition  number  of  the  closed') 

disp('     loop  acgf  matrix.  ') 

dispC    ') 

dispC    ') 

[acgf, bcgf ,ccgf ,dcgf ]=feedbk (algf,blgf,clgf,dlgf,2) 

tol=100*eps*norm( [acgf  bcgf]) 

[ acgf c, bcgf c,ccgfc,t,k]=ctrbf (acgf , bcgf , ccgf  ,  tol) 

tol=100*eps*norm( [acgf ; ccgf ] ) 

[acgf o, bcgf o, ccgf o, t ,k]=obsvf (acgf, bcgf, ccgf, tol) 

condacgf =cond (acgf) 

(strike  a  key  to  continue  ...)') 


of  the  closed  loop  system') 


dispC   ') 

dispC    ') 

dispC 

%pause 

%diary  hx29.dat 

dispC    ') 

dispC     Poles 

dispC    ') 

polcl=eig(acgf ) 

%hx29aly 

%hx29trsp 

end 
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format  short  e 

w=logspace (-3,3, 100) ; 

[magi, phase 1 ] =bode ( a lgf,blgf,clgf,dlgf, l,w) ;  magl=2  0*loglO (magi) 

[mag2 , phase2 ]=bode (algf , blgf , clgf ,dlgf , 2 ,w) ;  mag2=20*logl0 (mag2) 

[mag3 , phase3 ]=bode (algf , blgf , clgf , dlgf , 3,w) ;  mag3=20*loglO (mag3) 

semi logx(w, magi ( : , 1) ,w,mag2 ( : , 1) ,w,mag3 ( : , 1) ) 

title ('X-29  OPEN  LOOP  (GF)  BODE  PLOT  OF  "  alpha   "  ') 

xlabelC  FREQUENCY  -  rad/sec') 

ylabelCGAIN  -  dB') 

grid 

Imeta  h3bdl 

Ipause 

semilogx(w,magl ( : , 2) ,w,mag2 ( : , 2) , w,mag3 ( : , 2) ) 

title ('X-29  OPEN  LOOP  (GF)  BODE  PLOT   OF  "  q   "  ') 

xlabelC  FREQUENCY  -  rad/sec') 

ylabelCGAIN  -  dB') 

grid 

%meta  h3bd2 

%pause 

[magi , phase 1 ]=bode (acgf , bcgf ,ccgf , dcgf , l,w) ;  mag 1=20* log 10 (magi) 

[mag2 ,phase2 ]=bode (acgf, bcgf , ccgf, dcgf , 2 ,w) ;  mag2=2  0*loglO (mag2) 

[mag 3 , phase 3 ]=bode (acgf, bcgf, ccgf , dcgf , 3 ,w) ;  mag3=2  0*loglO (mag 3) 

serai logx (w,magl ( : , 1) ,w,mag2 ( : , 1) ,w,mag3 ( : , 1) ) 

title ('X-29  CLOSED  LOOP  BODE  PLOT  OF  "  alpha   "  ') 

XlabelC  FREQUENCY  -  rad/sec') 

ylabel( 'GAIN  -  dB') 

grid 

text  (15, -20,  ' TO  rl') 

text(  .22,-20,  ' TO  r2  ' ) 

text(.  1,-60,  ' TO  r3') 

meta  h3bd3 

pause 

semi logx (w, magi ( : , 2) , w,mag2 ( : , 2) ,w,mag3 ( : , 2) ) 

title ('X-29  CLOSED  LOOP  BODE  PLOT   OF  "  q   "  ') 

XlabelC  FREQUENCY  -  rad/sec') 

ylabel( 'GAIN  -  dB') 

grid 

text(  .1,-10,  ' TO  r2') 

text(.l,-30,  ' TO  rl') 

text  (.1, -58,  ' TO  r3') 

meta  h3bd4 
%pause 

dispC     ') 
end 
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APPENDIX  B 


State  Space  Realization  of  the  Open  Loop  Plant  G 


ag=[-.41 

.4427d- 

.547 

-.1164d 

0.0 


81d+00 
04   -. 
4d  +  01 
-01   - 
0.0 


0.0  0.0 
0.0  0.0 
0.0  0.0 
0.0  0.0 
0.0  0.0 
0.0  0.0 
0.0  0.0 
0.0  0.0 
0.0   0.0 

0.0   -.2697d 
0.0   0.0 

0.0   0.0   -. 

0.0   0.0 

-.1790d+02   0 


.9960d 
12d-04 
.3424d 
397d-0 
0   0.0 


-00 


145 


1479d+ 
3  0.0 
0  -.1 
Old+03 
0  0.0 
0.0 


+  00 
2   - 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 
04 

0. 
491d 

0. 

-.  6 


.  2269 
ld-05 
.2585 
509d 


d-01   -.1213d+00   -.1948d-01 
.2162d-04   -.3540d-05   0.0 


d  +  01 
03 

0 

1 

0 

0 

0 

0 

0 

0 

0 
0 


02d+05 
3d+03] 


bg=| 


.1479d+04  0.0; 
0.0  .1491d+04; 
.5308d+05   0.0] 


-.  1386d  +  01 
5362d-03 
0.0   0.0 


-. 1058d+01 
2912d-03   0 
0   0.0   0 
0 
0 
0 
0 
,0 
,0 
,0 
,0 


-.9493d- 

0.0   0. 

.3898d-0 

0  0.0  0. 

0.0 

0 

0 

0 

0 

0 


-.1143d+03   0.0   0.0   -, 


0. 

0. 

0. 

0. 

0. 

0. 
.Oe+04  0.0 
.0  1. Oe+04 
.0  0.0  1. 
2529d+01   0 


03 

0; 

2 
0 

0 
0 
0 
0 
0 
0 

0. 

0, 

Oe+04 ; 
.0  ... 


0; 
0; 


0.0   0.0   -.1149d+03   0.0   0.0   -.2536d+01 


0.0   0.0   -.1816d+04   0.0   0.0 


cg=[1.0 
0.0 


0.  0 
1.0 


0.0 
0.0 


0.0 
0.0 


0.0 
0.0 


0.0 
0.0 


0.0 
0.0 


0.  0 
0.0 


0.0 
0.  0 


0.0 
0.  0 


0.0 
0.0 


0.0 
0.  0 


0.0 
0.0 


0.0; 

0.0] 


dg=zeros (2 ) 
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APPENDIX  C 


State  Space  Realization  of  Reduced-Order  Controller  for  2  I/O  System 


acpr  = 

Columns  1  through  6 


■1 

.0399e- 

-02 

-1 

.22626- 

-05 

-1 

.04606-01 

4 

.80946-02 

-7, 

.68756-02 

1. 

,32746-02 

•1 

.12716- 

-04 

-1 

.01566- 

-02 

-4 

.67956-03 

8, 

.32526-02 

-3, 

.  18616-02 

-2, 

,80006-02 

•1 

.03406- 

-01 

-2 

.99776- 

-02 

-1, 

.95256+01 

2, 

.63946+01 

-3, 

.12806+01 

4, 

,92716+00 

•5 

,76156- 

-02 

6 

.01046- 

-02 

-2 

.70296+01 

-2. 

.56746+01 

4. 

,03966+01 

-7. 

,70336-02 

7, 

.86016- 

-02 

1 

.34196- 

-02 

3 

.17466+01 

3, 

.64836+01 

-1, 

,20306+02 

5, 

.35096+01 

•3, 

.40486- 

-03 

-3, 

.61386- 

-02 

-2. 

.74926+00 

5, 

.99976+00 

-2, 

,95726+01 

-9, 

.89846+00 

■4. 

.98646- 

-03 

-3. 

.97356- 

-02 

-3  , 

. 13026+00 

3. 

. 57646+00 

5. 

.92536+01 

-2, 

.47776+01 

3. 

.65356- 

-03 

-2, 

.91126- 

-02 

2, 

.90376-01 

1. 

.05136+01 

-2, 

.47176+01 

-2, 

.16856+01 

4. 

,79086- 

-02 

-1. 

.13006- 

-02 

1. 

.70106+01 

3. 

.28136+01 

-1. 

,47666+02 

1. 

.02906+01 

9, 

.20376- 

-03 

-1, 

.46526- 

-03 

3. 

.28776+00 

6. 

. 37586+00 

-2, 

,86616+01 

8. 

,90496-01 

8. 

.72006- 

-04 

6, 

.88446- 

-03 

-9, 

.96216-02 

-2, 

. 13296+00 

1, 

,64466+00 

3, 

,58466+00 

8. 

.30466- 

■04 

-3, 

.79216- 

■03 

1. 

.81966-01 

1. 

.40026+00 

-1, 

,96636+00 

-1. 

,89696+00 

8. 

,24456- 

-04 

1. 

.09426- 

-04 

3. 

.03026-01 

5, 

.13686-01 

-2, 

,52186+00 

1. 

,78646-01 

1. 

,77636- 

■04 

9. 

.00896- 

■04 

-3, 

.60946-02 

-3. 

.  19816-01 

4, 

,06686-01 

4. 

,54026-01 

2. 

31886- 

■05 

3. 

.65276- 

■05 

9. 

.57886-03 

6, 

.88326-03 

-7, 

,58966-02 

2. 

,28426-02 

Columns  7  through  12 


2. 

,58706-03 

9, 

.12866-03 

8. 

.38866-03 

8, 

.69406-03 

5. 

,30956-04 

1. 

,75036-03 

'4. 

,03366-02 

2. 

.52386-02 

-8, 

,40336-02 

-5, 

,50856-03 

-7, 

,00346-03 

-2, 

,65316-03 

9, 

,21186-01 

3. 

.49656+00 

2, 

,31616+00 

3. 

,  15646  +  00 

1, 

,54636-01 

6, 

,28546-01 

1. 

,01426+01 

-1. 

.37636+01 

1, 

,92466+01 

-3, 

.07666+00 

1, 

,85106+00 

-1, 

,36426-01 

5, 

,17796+01 

4. 

.64676+01 

1. 

,98536+01 

2. 

.47446+01 

-1. 

,16566+00 

4, 

.17216+00 

3, 

,22986+00 

1. 

.47706+01 

-4, 

,23326+01 

-7, 

,64186+00 

-3, 

,28326+00 

-2. 

.17056+00 

1, 

,39986+01 

3, 

.52546+01 

-6, 

,37786+01 

-5, 

,06286+00 

-4, 

,92786+00 

-2. 

.05586+00 

3. 

,70126+01 

-1. 

.43966+01 

9, 

.76796+01 

-3, 

,02196+00 

5, 

,00516+00 

5. 

,65936-01 

2, 

,04926+01 

-7. 

.64986+01 

-2, 

,43056+02 

-2. 

,94396+01 

-4, 

,12336+01 

-1. 

,79446+01 

1. 

,56486+00 

-4, 

.63696+00 

-8. 

,70736+01 

-4. 

,56006+01 

-2. 

,46196+01 

-2. 

,10836+01 

4. 

,47926+00 

6. 

.67106+00 

1, 

,68076+01 

2, 

,01426+01 

-2. 

,72376+01 

-2. 

,58956+01 

2, 

,33566+00 

-3, 

.89246+00 

-1, 

,12756+01 

-9, 

,79366+00 

3. 

,29616+01 

-1, 

,25506+02 

2. 

,79696-01 

-2. 

.15206-01 

-8, 

,10786+00 

-8, 

,05346+00 

1. 

,93346+00 

-2. 

,75686+01 

5, 

,60356-01 

9, 

.16706-01 

2. 

, 50246+00 

2, 

,23966+00 

-7, 

,30706+00 

5. 

,73806+01 

3, 

.03106-02 

2. 

.54826-02 

-2, 

,01056-01 

-2. 

,07316-01 

-2, 

,03026-01 

1. 

,17006+00 

Columns  13  through  15 


•4.87876-04 
-1.17986-03 
•1.86416-01 

6.21646-01 
•1.90866+00 
•1.92156-01 
•7.37426-01 

1.40356+00 
-5.73136+00 

2.34046+00 
•8.15436+00 

2.46476+00 
■4.56316+01 

3.33966+01 
■2.70146+00 


-4.76006-04 
-4.62236-04 
-1.77856-01 

3.97056-01 
-1.59496+00 

9.88206-02 
-2.48866-01 

8.09526-01 
-1.74346+00 

3.53056+00 
-2.59546+00 

1.51066+01 
-6.34406+01 
-7. 13156+01 
-3.07936+00 


•1. 16486-05 
-5.06366-05 
-4.58156-03 

2.16766-02 
■5.43296-02 
-1.39226-02 
■3.29506-02 

5.17506-02 
■2.63086-01 

1.14766-02 
■3.65116-01 
-1.07586+00 
■3.31306+00 
•6.29576+00 
-1.44186+02 
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bcpr  = 


5. 

25026+00 

6. 

46186-01 

3. 

28306-01 

2. 

,91016+00 

2. 

59706+01 

7. 

,30916+00 

1. 

.55898+01 

-7. 

,12556+00 

1. 

,96096+01 

-4. 

,  14076  +  00 

2. 

,  18936-01 

5. 

,28166+00 

5. 

,57836-01 

5. 

,84526+00 

■1. 

.43666+00 

4. 

,05636+00 

1. 

.23106+01 

2. 

.59796-01 

•2, 

.35226+00 

-5. 

.26536-02 

3, 

.41656-01 

-9, 

.60276-01 

-2 

.76706-01 

5 

. 18996-01 

-2 

,06466-01 

-3 

.91676-02 

6 

.07656-02 

-1 

.23866-01 

-5 

.21776-03 

-5 

.88836-03 

ccpr  = 

Columns  1  through  6 

-3.87796+00    1.98296+00   -1.92576+01    7.12366-01   -1.1251e+01    5.2413e+00 
-3.59796+00   -2.15516+00   -1.8895e+01    1.71266+01   -1.65856+01   -6.86776-01 

Columns  7  through  12 

4.45016+00   -7.61606-01    9.8134e+00    2.16756+00    7.86126-01    5.8776e-01 
-3.83076+00    4.23526+00   -7.43606+00    9.15146-01   -6.48746-01    2.10196-02 

Columns  13  through  15 

2.43806-02   -4.37676-02    2.78526-03 
-2.08726-01   -1.30836-01   -7.35806-03 


dcpr  = 

-3.98226-01   -9.29516-02 
-6.39886-02   -1.49366-02 
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APPENDIX  D 


Closed  Loop  State  Space  Realization  of  Controller/Plant  Series 


Columns  1  through  6 


-1.0399e-02 
-1.1271e-04 
-1.0340e-01 
-5.7615e-02 

7.8601e-02 

-3.4048e-03 

-4.9864e-03 

3.6535e-03 

4.7908e-02 

9.2037e-03 

-8.7200e-04 

8.3046e-04 

8.2445e-04 

-1.7763e-04 

2.3188e-05 

■1.2486e+00 

6.2283e-01 

■2.8230e+00 

4.5916e-01 

■2.6795e+00 

3.3498e+00 

-1.1953e+00 

-1.9331e-01 

1.1620e+00 

-9.4959e-01 

■4.6009e-01 

1.4474e-01 

1.3788e-02 

1.3630e-02 


-1.2262e-05 

-1.0156e-02 

-2.9977e-02 

6.0104e-02 

1.3419e-02 

-3.6138e-02 

-3.9735e-02 

-2.9112e-02 

-1.1300e-02 

-1.4652e-03 

6.8844e-03 

-3.7921e-03 

1.0942e-04 

9.0089e-04 

3.6527e-05 

-2.9277e-01 

-3.5136e+00 

4. 3606e+00 

-1.8217e+00 

7.0049e-01 

-6.0600e-01 

1.4172e+00 

9.0043e-01 

-1.4093e-01 

5.5372e-01 

1.0922e-02 

7.9646e-02 

-9.7163e-02 

-6.9743e-03 


-1.0460e-01 
-4  .67956-03 
•1.9525e+01 
■2.7029e+01 

3.1746e+01 

■2.7492e+00 

-3. 1302e+00 

2.9037e-01 

1.7010e+01 

3.2877e+00 

-9.96216-02 

1.81966-01 

3.03026-01 

•3.60946-02 

9.57886-03 

■6.44006+00 

2.27046+00 

-1.32686+01 

1.87166+00 

■1. 34786+01 

1.69206+01 

-5.72856+00 

-7.53606-01 

5.88706+00 

-4.69816+00 

-2.34256+00 

7. 58346-01 

4.52746-02 

6.76866-02 


4.80946-02 

8. 32526-02 

2.63946+01 

■2.56746+01 

3.64836+01 

5.99976+00 

3.57646+00 

1.05136+01 

3.28136+01 

6. 37586+00 

-2.13296+00 

1.40026+00 

5.13686-01 

-3. 19816-01 

6.88326-03 

4 .06756+00 

1.30556+01 

■1. 15046+01 

6.45616+00 

3.25236+00 

■5.17766+00 

■3.10236+00 

■3.26836+00 

-2.08166+00 

■1.06416-01 

1.00926+00 

■6. 59916-01 

3.68876-01 

-2.48596-03 


-7.68756-02 

-3.18616-02 

-3. 12806+01 

4.03966+01 

-1.20306+02 
-2.95726+01 

5.92536+01 
-2.4717e+01 
-1.4766e+02 
-2.8661e+01 

1.64468+00 
-1.96636+00 
-2.52186+00 

4.06686-01 
-7.5896e-02 
-5.05536+00 
-3.  10836+00 
-3.70326+00 
-1. 10916+00 
-8.80466+00 

1. 14226+01 
-2.22836+00 

6.72296-01 

4 .06876+00 
-2.65046+00 
-1.68016+00 

6.56366-01 
-9.86236-02 

3.95416-02 


1. 32746-02 
-2.80006-02 

4.92716+00 
■7.70336-02 

5.35096+01 

-9.89846+00 

-2.47776+01 

-2.16856+01 

1.02906+01 

8.90496-01 

3.58466+00 

-1.89696+00 

1.78646-01 

4.54026-01 

2.28426-02 

3.93916-01 

-5.28056+00 

7.86796+00 

-2.82516+00 

2.69136+00 

-2.98996+00 

2.73536+00 

1.37486+00 

-9.40876-01 

1.37816+00 

3.10206-01 

1.78306-02 

-1.43826-01 

-1.84296-02 


Columns  7  through  12 


2, 

,58706-03 

9, 

,12866-03 

8. 

,38866-03 

8, 

,69406-03 

5. 

,30956-04 

1. 

,75036-03 

4. 

,03366-02 

2, 

,52386-02 

-8, 

,40336-02 

-5, 

,50856-03 

-7, 

,00346-03 

-2. 

,65316-03 

9, 

,21186-01 

3. 

,49656+00 

2. 

,31616+00 

3, 

,15646+00 

1. 

,54636-01 

6. 

,28546-01 

1, 

,01426+01 

-1. 

,37636+01 

1, 

,92466+01 

-3, 

,07666+00 

1, 

,85106+00 

-1, 

,36426-01 

5. 

,17796+01 

4. 

,64676+01 

1, 

,98536+01 

2, 

,47446+01 

-1, 

,16566+00 

4. 

,17216+00 

3, 

,22986+00 

1. 

,47706+01 

-4. 

,23326+01 

-7. 

,64186+00 

-3. 

,28326+00 

-2. 

,17056+00 

1. 

,39986+01 

3, 

,52546+01 

-6. 

,37786+01 

-5. 

,06286+00 

-4. 

,92786+00 

-2, 

,05586+00 

3, 

,70126+01 

-1. 

,43966+01 

9, 

,76796+01 

-3, 

,02196+00 

5, 

,00516+00 

5, 

,65936-01 

2, 

,04926+01 

-7, 

,64986+01 

-2, 

,43056+02 

-2, 

,94396+01 

-4, 

,12336+01 

-1. 

,79446+01 

1, 

.56486+00 

-4. 

,63696+00 

-8, 

,70736+01 

-4. 

,56006+01 

-2, 

,46196+01 

-2. 

,10836+01 

4, 

.47926+00 

6, 

,67106+00 

1. 

,68076+01 

2, 

,01426+01 

-2, 

,72376+01 

-2. 

,58956+01 

2. 

.33566+00 

-3, 

,89246+00 

-1. 

,12756+01 

-9, 

,79366+00 

3. 

,29616+01 

-1. 

,25506+02 

2. 

.79696-01 

-2, 

,15206-01 

-8, 

,10786+00 

-8. 

,05346+00 

1, 

,93346+00 

-2. 

,75686+01 

5, 

.60356-01 

9, 

,16706-01 

2, 

,50246+00 

2, 

,23966+00 

-7, 

,30706+00 

5. 

,73806+01 

3, 

.03106-02 

2, 

,54826-02 

-2. 

,01056-01 

-2, 

,07316-01 

-2. 

,03026-01 

1, 

,17006+00 

•4. 

.22606-01 

9, 

,06776-01 

-6. 

,96116-01 

4. 

,42436-01 

-6, 

.81366-02 

6. 

,70256-02 

■7, 

.08096+00 

4. 

.07496+00 

-1. 

,48066+01 

-1. 

,22476+00 

-1. 

,22856+00 

-5. 

,13756-01 

9. 

,05166+00 

-4. 

, 16306+00 

1, 

,92226+01 

2, 

,37826+00 

1. 

,57866+00 

8. 

,10736-01 

•3. 

.68876+00 

2, 

,05336+00 

-7, 

,73266+00 

-6, 

,91986-01 

-6. 

,40516-01 

-2. 

,77876-01 

1. 

.74066+00 

3. 

,02176-01 

4. 

,00806+00 

1. 

,31406+00 

3. 

,12176-01 

3. 

,18246-01 

1, 

.63876+00 

-7, 

,11456-01 

-3, 

,89386+00 

-1, 

,56876+00 

-2, 

,97216-01 

-3, 

,62456-01 

2 

.97766+00 

-1, 

,23186+00 

6, 

,36226+00 

8. 

,89256-01 

5, 

,20366-01 

2, 

,86966-01 

1. 

.81906+00 

-1, 

,02966+00 

3. 

,80826+00 

3, 

,27956-01 

3, 

.15716-01 

1, 

,34506-01 

•4, 

.30376-01 

-3, 

,32506-01 

-1, 

,06386+00 

-5, 

,25156-01 

-7, 

,91976-02 

-1. 

,16636-01 

1. 

.22556+00 

-2, 

,70796-01 

2, 

,68526+00 

5, 

,49476-01 

2. 

,16016-01 

1, 

,52876-01 

8. 

,09926-02 

1, 

,8717e-01 

2, 

,35416-01 

1. 

,95636-01 

1. 

.58776-02 

4. 

,02916-02 

1. 

.40066-01 

-1, 

,61666-01 

2. 

,69956-01 

-3, 

,87516-02 

2. 

,36666-02 

-1, 

,77136-03 

•1, 

.95376-01 

1, 

,14186-01 

-4, 

,08016-01 

-3, 

,24286-02 

-3. 

,38816-02 

-1. 

,39176-02 

•1, 

.56516-02 

2. 

.68246-03 

-3. 

,45126-02 

-7. 

,62006-03 

-2, 

,76476-03 

-2, 

,06656-03 
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Columns  13  through  18 


4 

8787e-04 

-4 

76006-04 

1 

1798e-03 

-4 

.62236-04 

1 

8641e-01 

-1 

77856-01 

6 

.21646-01 

3 

.97056-01 

1 

9086e+00 

-1 

.59496+00 

1 

9215e-01 

9 

.88206-02 
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1 
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-1 
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3 
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-5 
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7 
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8 
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7 
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9 

19806+00 

5 
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6 
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1 
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-1 
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7 
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8 

36576-01 

2 

87186-02 

-4 

92666-02 

1 

73196-01 

-1 

99336-01 

7 

29816-03 

-8 

98186-03 

2 

23086+00 

8 

29046-02 

6 

61926-01 

-1 

71406+00 

1 

53086+00 

-1 

44316+01 

5 

66146-01 

1 

20236+01 

2 

07136-01 

-7 

73576+00 

2 

07576-01 

7 

78286+00 

4 

51336-01 

-9 

86846+00 

2 

79196-01 

-5 

45536+00 

5 

07256-02 

2 

19146+00 

1 

80426-01 

-4 

34996+00 

5 

82826-03 

-6 

09376-01 

2 

27956-02 

-2 

96106-01 

3 

04016-02 

5 
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2 

21336-03 

5 

67986-02 

Columns  19  through  24 


-3.64266-01 

1.82696+00 

7.13066-01 

-6.61206+00 

3.04226-01 

3.19796+00 

3.49366+00 

2.68816+00 

1.93936+00 

3.08036-01 

-6.36606-01 

3.57376-01 

5.91276-03 

-8.45246-02 

-2.84596-03 
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-1. 19606+01 

-3.27526+01 

1.24766+01 

-9.81556+00 

4.19196+01 

3.62046+01 

-8.93476+00 

2.27626+01 
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-2.15956+00 

■1.13596+00 

■1.24536+01 

■2.51626+00 

8.80936+00 

-2.35356+00 

■2.71786+00 

■1.23266+00 

4.29866+00 

8.65466-01 

2.91246-01 

■1.24426-01 

9.08446-02 

3.15816-02 

4.40596-03 

■1.77556-01 

2.85406-01 

8.16716+00 

2.54036+01 

-3.36646+01 

6.67746+01 

■3.79026+01 

-2.33856+01 

2.72186+01 

■3.04656+01 

-9.70246+00 

1.60346+00 

1.99876+00 


2.46836-01 

1.94996+00 

3.96056+00 

-5.29156+00 

-2.07256+00 

3.51046+00 

3.87406+00 

2.75226+00 

6.03526-01 

4.71846-02 

-6.51596-01 

3.55386-01 

-1.88696-02 

•8.46266-02 

-3.73986-03 

1.85916-02 

-6.78166-01 

-1.22786+01 

-6.52446+01 

3.11726+01 

-1.14996+02 

-6.78436+01 

9.09796+01 

■8.51956+01 

8.63836+01 

2.93936+01 

-6.09566+00 
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4.05756-01   -1.10276+00 
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Columns  2  5  through  29 


7.6382e-01 

6.4187e-01 

4.7395e+00 

1.5408e-01 
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1.2600e+00 

1.4313e+00 

7.8041e-01 
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1.61186+00 
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3.20246+00 

9.22556+00 

-9.29926-01 
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bcgf  = 


5.25026+00 
-3.28306-01 

2.59706+01 

1.55896+01 
-1.96096+01 

2.18936-01 

5.57836-01 
-1.43666+00 
-1.23106+01 
-2.35226+00 

3.41656-01 
-2.76706-01 
-2.06466-01 

6.07656-02 
-5.21776-03 
-5.70086-02 

3.08296-01 
-5.12966-01 

1.68506-01 
-2.23956-01 

2.59356-01 
-1.84386-01 
-8.11486-02 

8.46666-02 
-1.02736-01 
-3.00926-02 

3.11366-03 

8.30686-03 

1.40016-03 


6.46186-01 

2.91016+00 

7.30916+00 

-7.12556+00 

-4.14076+00 

5.28166+00 

5.84526+00 

4.05636+00 

2.59796-01 

-5.26536-02 

-9.60276-01 

5.18996-01 

-3.91676-02 

-1.23866-01 

-5.88836-03 

-1.33066-02 

7.19606-02 

-1.19736-01 

3.93316-02 

-5.22746-02 

6.05366-02 

-4.30386-02 

-1.89416-02 

1.97626-02 

-2.39786-02 

-7.02396-03 

7.26776-04 

1.93896-03 

3.26796-04 
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ccgf  = 

Columns  1  through  6 


Columns  7  through  12 


Columns  13  through  18 


Columns  19  through  2  4 


1.4464e-01 
-6.1148e-01 


3.5829e-01    3.4970e-02 
4.3076e-01   -6.6611e-01 


Columns  2  5  through  2  9 

-1.1672e-01    4.0736e-02 
-2.3373e-01   -7.6955e-02 


0   -6.4236e-01 
0   -8.3474e-02 


1.1544e-01 
-1.3417e+00 


6.0999e-02 

-1.5794e+00 


-2.1996e-02    7.9072e-02   -2.7843e-02 
5.3449e-01   -2.9948e-01    2.2334e-01 


3.0755e-02   -3.1054e-02   -1.7417e-03 
2.3094e-02   -1.1998e-02    3.0561e-03 


dcgf  = 

0 
0 


74 


APPENDIX  E 


Poles  of  The  3  I/O  Reduced  Order  Closed  Loop  System 


-9.1468e+03 

-5.5159e+02  +  5.3384e+02i 

-5.5159e+02  -  5.3384e+02i 

-2.2777e+02  +  2.3189e+02i 

-2.2777e+02  -  2.3189e+02i 

-1.4457e+02 

-1.4319e+02 

-1.1949e+02 

-8.9079e+01  +  1.8937e+01i 

-8.9079e+01  -  1.8937e+0li 

-1.7935e+01  +  6.0318e+01i 

-1.7935e+01  -  6.0318e+01i 

-5.2531e+01  +  4.8268e+01i 

-5.2531e+01  -  4.8268e+01i 

-4.9021e+01  +  4.9207e+01i 

-4.9021e+01  -  4.9207e+01i 

-3.1654e+01  +  9.6963e+00i 

-3.1654e+01  -  9.6963e+00i 

-2.9006e+01 

-2.0127e+01 

-8.7118e+00  +  9.5608e+00i 

-8.7118e+00  -  9.5608e+00i 

-8.0945e+00  +  6.1691e+00i 

-8.0945e+00  -  6.1691e+00i 

-1.4512e+01 

-1.8891e+00 

-1.9503e+00 

-1.0702e+01 

-2.7168e+00 

-8.6993e+00 
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